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INTRODUCTIOH 

ASM80 is a tv;o pass symbolic assembly program for the IJ3TSL 
8080 Microprocessor. The assembler runs on the NOVA line of mini- 
computers under DOS or RDOS. 

ASM80 accepts symbolic assembly language statements as input, 
producing absolute binary as output, . The output is suitable for 
input to the FROM 1702A programmer, and the standard binary loader 
for the 8080 system. 

In describing statements in this manual, the follov;ing convsn- 
. 1 

tions vrill be used: 






1) All numbers 
wise . A number 
that the number 
is 16 (decimal) . 
ing letter K as 


are decimal unless explicitly marked other- 
follo'-,;ed by the characcor X or H, indicares 
is expressed in the hexidecimal radix, as Ifx 
Octal numbers are indicated by a terminar- 
the final character as Ip'K, which is 8 decimal. 


2) In statement descriptions, items enclosed in square brackets 
are optional . 

3) Items enclosed in brackets "{}" are optional and 

may be repeated any number of times. 


4) The symbol CR indicates a carriage return character. 


5) In those cases where spacing characters are mandatory 
in a statement, the symbol *'A" will be used to indicate any 
number of blanks or tabs in the line. Unless so noted, spaces 
nre not significant. 


1,1 ASM80 Language 

ASM80 prograras are generally prepared on a NOVA system, 
using the standard editor. Any other source of ASCII input, 
which can be loaded on a NOVA system as a text file, and v:hich 
conforms to the programming conventions outlined in the section 
is acceptable as input to the Assembler. 

Each statement is written on a single line, terminated 
by a carriage return. Blank lines are allowed so the user 
may group statemshts. Statements are format free; that is, 
individual, elements of a statement need not appear in parti- 
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cular columns of the input line. The elements of statements 
are defined by specific delimeter characters. 

There are four tvpes of elements in an ASM 80 statement. 
The assembler identifies the element type by its relative 
location, and by the presence of special delimeters. The in- 
dividual elements need not appear in specific columns, but 
their order is fi>:ed. Statements have the general form: 


label ; operator opsrandl ,operand2 ; comments (carriage return] 


Each statement is proceeded by the assembler as required for 
given by the programmer. Some operators result 
in the generation of binary mnchina code, v;hile others are^ 


simply interpreted by the assembler to control the assembly 
process . 


The second type of ASMSO statement does not have an 
operator, but rather consists of a list of expressions. In 
general, the format of such an expression statement is: 


label t expression! r expression2 express icnN ; comment (carriage 

return) 


. Expressions will be discussed in detail later, but they 
may be as simple as literal numbers. . An example of this would 

be: 

BITtlSK: 01X,02X,04X,08X,10X,20X,40X,80X 


This statement generates 8 bytes of data, corresponding in 
sequence to the 8 simple expressions in the above statement. 
The symbol 3ITMSK is assigned the value of the byte address 
of the first byte in this table. 

By default, the assembler generates a single byte of out- 
put for each expression in the list. The psuedd-op "EXPISK” 
allows the. programmer to change the default output length to 


ORli&lNAIi PAGE IS . 
OF POOR QUAIlHi 



1,2 (3 r or 4 bytes. See Section 


for details. 

The third type of statep.en-t allows the user to assign 
a value to a symbol without using the symbol as a label. This 
statement is the assignment statement, and the general form 
is: 

syroboll , syp.bol2 , . . . symbplN « expression ? comment (carriage return) 

The assembler simply evaluates the single expression given, 
then assigns the value of the expression to the symbols given 
in the Isit to the left of the equals sign. An example of 
a possible expression is: 

TOP . USED . COFJB =* TOP , USE D . COBE + 8 

This statement reassigns a nsv; value to the symbol "TOP. USED. 
CORE", which is its old value plus 8, 

The symobls v/hich are assigned a value may not be used 
as labels in the program since this would provide conflicting 
Information about the location of the statement specified by 
the label, 

1.1.1 Labels ' 

A label is a symbol, created by the programmer to 
identify a particular statement in his program. A label 
must be a properly formed symbol (see section 1.2) fol- 
lowed by a dolOn, VThile our examples of a statement 

have only showed one label, it is legal to define multi- 
ple labels for a single statement, 

label 1 : label2 ; label3 : operator ope rmiacan (carriage Re tui^> 

In practice, the assemblar assigns a va3^ to a label 
when it is defined. This value is the byte address, in 
INTEL memory, of the first by t^ of the binai^ generated ^ ^ 
by the statement following thef labei. " 


V 


is assigned a value between U and FFFFX iiiGlusive, 


1.1.2 Operator 


An operator may be one of the INTEL 8080 maehine 
instructions defined in the specifications for the 8080 
microprocessor, or a command for the assembler. A 
command for the assembler is often referred to as a 
psuedo-operation , or psuedo-op. 


The Operator determines the format of the remain- 
der of the statement, such as the number and type of oper- 
ands • 


1.1.3 Operands 


Operands are items which specify the data item, 
or the address of the data item, to be operated on by 
an operator. Operands may be as simple as numbers, or 
may be eKprcssions. Generally operands that refer to 
the address of a data item are symbolic. The symbol is 
the label assigned to the data item by the programmer. 

If an operator requires multiple operands, they are simply 
listed after the operator— separated by commas. 


Consider the following examples of actual 8060 statements; 


Statement 


le of Operand ( s ) 


ADI 16. 


S^«A NEWVALUE 


STA NEWALUE+1 


Literal Decimal Constant 
Symbolic Reference to a Label 
Symbolic Reference to the byte 


mi B, THRESHOLD/SCALE 


just 1 byte after the label ''NZKVAL'JI” 

The first operand is a symbolic ref e 1*7 
ence to 8080 register B (see sec tier, 
for predefined si-rntbpls) . The second 
operand is an expression, v/hbse val-a 
will become the iirmediate byte valvs 
for the move immediate instructicn. 






-5- 

ASM80 will accept an expression for any operand. The 
value of the expression will be tested to determine that 
it is in the legal range of valid operands for the opera- 

f 

ter it is used with* 


1.1.4 Coininents 

The progranmer may piece notes about the program 
on every statement line by marking the beginning of the 
coiTiinent v?ith e* semicolon. The assembler will ignore any 
text after a semicolon. The programmer may use an entire 

line as a comment, by starting the line with a semicolon . 

# 

To make it convenient to enter large blocks of 
comm.ants, a special type of Gcir;mcnt called the block 
comment is allox'v’ed. All the text contained in the special 
delime terpairs “/*" and "*/" v^ill be ignored by the asssiri 
An example of their use is : 

/* any text the prog ramir.er desires ..... . 


• • « -A • * 

*/ 


Block comments have the additional property that 
they nest, thus the assembler will allow a block comment 
to contain another block comment. This allows the pro- 
grammer to quickly block out a section of a program with^ 
out worrying about enclosed comments. For example; 


/♦ DELETi:iG INITIALIZATION FOR TESTING 
Statements 

/*this section of code initializes the slot control qUe 

1) Scanning for ..... 

2) Initializing the map table* . • . 


• • « 






stateTRents 


END 0? DEIiETICN FOR TESTING */ 

In this example, the assembler recognizes the enclosed 
block co.niment as such, and. does not terminate tiie enelos- 
iiig comment prematurely. 

1,2 SY^NS0LS 

The programmer may create symbols for use as labels or 

operands. A symbol consists of from one to 31 characters chosen 

from the following set: , 

* 

The 26 letters A,B,C, . . , ,vr,X,Y,Z 

the 10 digits 0,1,2,3, . . .Q,9 

3 special characters Dollar sign ($) , percent (%) , and 

period ( , ) 

The first character of a symbol may not be a digit. Examples 
of legal and illegal symbols are: 



Legal 

A 

A1 

%PSTATUS 
TOP. USED. CORE 



(Digit illegal as 1st char 
(No spaces^ in si'mbcls) 

(& not in legal set) 


If a symbol is longer than 31 characters, the excess characters 
are truncated, thus two symbols with the same initial 31 charac- 
ters are indistiguisable to the assembler. 
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1.2*1 Special Reserved Symbols 



A number of symbols have permanent values assigned 
to them by the assembler. These symbols include; 


1) All 8030 instruction syirliols, 

2) All ASM80 psuedo-oporators 


For convenience in referring to the 8080 CPU registers. 


the following symbols have been predefined to select the 
corresponding regi s ter , 


Symbol 

A 

B 

C 

D 

E 

H 

L 

M 

SP 



PSW 


Register 

7-CPU accumulator 
d-CPU register B 

1-GPU register C / 

2 

3 

4 

5-CPU register L ' 

6“Specifies K,L pair 

In certain instructions, selects the 
stack pointer 

In PUSH, POP instructiens , selects 
the A register and Program 
Status t'ford. 


In addition to the above symbols, the symbol is 
used to reference the current byte counted of the assembler. 
It may be used in an expression as any other symbol vrould 
be. At each ocourrence, its value is the value of the 
byte counter at the start of the statement. 


1.3 NUMBERS 


Literal numbers may be used in statement to represent 
specific values. The assembler interprets a number in the 
source based on the current input radix and any local radix 
specifiers that may modify the number. A simple integer num- 
ber. may consist of characters from the following set: 

« 

the digits 0|lf;2, ... ,8,9 
six letters A,B,G,b,E,P 


\ 


1 






8 


A number must begin with a digit. The simple integer num- 
• Trtsiif Krs TirtHS ■Fi a nridix snecif ier drawn from the 

•***-"“ •' -“'U ^ ^ 

following set: 


Specifier 

Kadix 


• 

Decimal 


K 

Octal (8) 


d 

Octal (8) 


B 

Hexidecimal 

(16) 

X 

Hexidecimal 

(16) 


Rn - 

Base n, where ri is in the 

range 

2 to 16 


inclusive, n is always 

given 

in the 


. decimal base . 



The specifier is appended to the number 1,3, 



Number 

Value in Decimal 



lOK 

8 



lOX 

16 



-10R2 

»-2 



OFX 

15 



By default, if no local 

radix specifier is given. 

the assembler 



assumes a number which is expressed in the decimal radix. The 
user may change the default by using the "IRDX” psuedo-operation. 
See Section 1.8.3 for details of this psuedo-op. In any ease, 
a local radix specifier always overrides the default input 
radix. ' 





1.4 STRING TOKENS 
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1.5 EXPRESSIONS 

An expression is a series of numbers and/or symbols 
separatisd. by a special class of delimeters .called arithmetic 
or logical operators. (We will' simply refer to arithmetic 
or logical operators as "operators”— there should be no con- 
fusion with operators that are 8Q80 instructions.) 'The 
assembler evaluates the expression at assembly time by ex- 
plicitly carrying out the operations indicated on the values 
of the numbers and sj'mbols. The following operations are 
accepted by the assembler: 


Opgy^^'hions 
NOT 
* , / 

+ r - 
& 

1 r XDR 


Meaning 

logical complement 
multiplication. Division. , 
addition, subtraction 
logical and 

logiGal or, logical exclusive 


All arithmetic operations are signed integer 32 bit operations. 
W1 logical operations are also carried out over 32 bit values* 
This is possible since the assembler al^cays maintains a 32 bit 
value for every symbol or number, even if that precision is 
not required for its particular occurrence. For example, 
labels are kept as 32 bit numbers, but in practice the value 
of a label must always fit in 16 bits since the 8080 can only 

handle memory addresses of 16 bits. At any point where the 

* 

actual value must be truncated to fit the application^ as an 
address field, the assembler checks to verify that the value 
is in the legal range for that application. 

The order in which the assembler evaluates an expression 
is determined by the relative precedence of the operators in 
the expression. In the previous list of operators, the opera- 
tors were listed in descending precedence , with those operators 
of equal precedence on the same line. Thus in evaluating the 
expression: X = A+3’^C, the assembler will perform the multi- 

0 Poor 


N 
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pliciition first, then add A to the product of 3 and C. 

If operators of equal precedence occur in the sanie expression, 
they are evaluated left to right. The following table of 
expressions derionstrates precedence evaluation. 


Expression Value 

5+3/2 6 

5+2/3 5 

5*2/3 3 

5/^2*3 6 

OFFX&l*2 2 


For convenience, the assembler- allows the use of paren- 
theses to override the normal precedence evaluation of expres- 
sions. Any expression enclosed in parenthesis will be evalu- 
ated before applying the operation (s) just ^outside the paren- 
thesis. For example, in the expression: (A+i3)*C, the asser.blar 
will perforin the addition first, since it must evaluate any 
enclosed expression before it "knov.’s" the value to use in the 
multiplication operation , 

1.6 ADDRESS ASSIGNI-ENTS 

As source statements are processed, instruction and data 
bytes are assigned to consecutive memory addresses. The byte 
counter is incremented as required by a statement to assign 
its data bytes to memory. Some statements, such as assign- 
ment statements do not generate data and thus do not affect 
the byte counter. 

1.6.1 Referencing the Byte Counter 


The special symbol ' . ' (simple period) may be used 
in any expression to reference the value of the byte counter. 
The value of the byte counter used is the value for the 
first byte of the data generated by the statement, regard- 
less of the number of bytes generated. Thus the symbol 
*•' has the value that would be assigned to a label, if 
one occurred on that line. For example. 


**LAEEL: JMP LABEL" and 

are equivalent instructions. 


"JtiP." 








12 


1.7 INSTRUCTION FOPi>!ATS' 




OiRKaKAi; PA6£! 13 
op POOS QO^UXS 
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1.8 ASM80 PSUE DO-OPERATORS 

A SDGcial class of oredefined psuedo-operators are used 
to control the operation of the assembler. The general for- 
mat of these psuedo-operations is the same as for true opera- 
tors • 

1.8.1 END Statement 

The END statement marks the end of the program 
for the assembler. Any program statements after an 
END statement are ignored by »,SM80, An optional operand 
allows the programmer to specify the starting address 
of the program. For . example; 

END expression 

The assembler evaluates the expression, and if 
binary output is veing produced, an end block is generated 
with the expression value as the starting address of the 
program. The loader for the 8080 system will use this 
address to start the program executing. If the expression 
is missing, a value of -1 (FFFFX) is output for the start- 
ing address. By convention a -1 will stop the loader, 
but not start the program executing. 

1.8.2 LOG Statement 

The LOG statement allov/s the programmer 'to force 
the program byte counter to a new value. This allows 
altering the normal assignment of statements to the 
next sequential location. The LOG statement takes the 
form; 

label ; LOG expression ; comment (Carriage return) 

V 

The expression is evaluated, the byte counter is set to 
the value. If the byte count is decreased as a result 
of the LOG statement a v/arning message is output, but 
it is accepted by the assembler. The value of any label 
on an LOG statement is the byte address that would have 
been assigned to an instruction statement or expression 
statement if one had been there. Thus the LOG statement 


1 
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may be used to allocate tables, ivith the label taking on the 


1.8.3 IRDX Statement 

The IRDX statement allows the programmer to change 
the default iiiput radix for conversion of numbers. The 
general form of the IRDX statement is: 

IRDX expression ; comment (carriage return) 

The expression is evaluated, and if the value is in the 
range 2 to 16 inlcusive then the default input radix is 
set to the value. Not that the expression will be evalu- 
ated in the current input radix (before the IRDX is com^- 
pletely processed) so -che user must use caution in speci- 
fying numbers in the expression. It is a good practice 
to always use a local radix specifier ' to assure you get 
the radix you desire. For example IRDX 8." will alv/ays 
set the input radix to octal regardless of the previous 
input radix . 

1.8.4 ORDX Statement 

The ORDX statement allows the programmer to change 
the default output radix of the assembler. Normally the 
assembler outputs program listings with the insvtruction 
locations and values in hexidecimal. The general form 
of the ORDX statement is: 

ORDX expression ; comment (carriage return) 

The expression is evaluated, and if the value is 2,8,10, 
or 16, then the output radix is changed to this new value. 
The listing v;ill be formulated as required for the speci- 
fied output radix. 



2.0 ASM0O OPERATION 


ASM80 is designt 2 d to execute on a NOVA line computer system, 
under DOS or RDOS. It is assumed that the reader is generally 


Jl^£UllXXXOir WXuH w*iw 


2.1 Creating the Source Program 


Source program for the ASM&0 can be created using the 
standard KOVA editor. As stated in earlier sections, each 
statement must be contained on a single line of source. A 
line may have as many as 132 characters, terminated by a earriage 
return • 

It is often convenient to break a program up into pages 
organized by the logical functino of the program. ASM30 
accepts multipage programs as input, and if a listing is de- 
sired each page of input will be started on a new pace of the 
listing file. 

2.2 Executing the Assembler 

ASK80 is invoked under DOS /RDOS as any other program 
would be; that is, the user simply types in toh name, "AS’lSO'' 
follov;ed by the commands to the assembler. ASM80 will process 
the corcnand line, then perform the indicated functions. 

2.2.1 Basic ASM80 Command Line ^ 

The basic command, under DOS/RDOS, to invoke the 

assembler is: 



ASM80 filename (carriage return) 

This command will assemble the program filename . 3y 
default# no listing of the program will be produced. 

If any error occurs, the error messages will appear on 
the system console device ($TTO) . An absolute binary 
output file, in hexldecimal format, will be produced in 
the file "f ilename .Bg'*. Any previous version of " filename .BS" 
will be deleted. 


If a programmer desired, his program may be broken 
up into several files. In this ease he should list the 
the filenames on the command line in the order they are 
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to be. read by ASM80. As an example consider assembling 
a 3 file program consisting of a parameter file, PAR!I: 
a standard system initialization file INIT; and the users 
program, SCAijNER, The command line for this operation 
would be ; 

ASM80 FARM INIT SCJ^NEH (carriage return) 

This command will correctly assemble the program, but 
the binary output will go into a file called "BARI*i.B8", 

It is not likely that the programmer actually desires 
this, so the assembler accepts command line switches 
which modify its normal operation for the users con- 
venience. Switches are discussed in the next, section^ 
but briefly the following command line would cause the 
binary output to go the file "SCAMMER. B8" for the user. 

ASM80 PARl-I INIT SCAI^NER SCAISMER. B3/B (carriage return) 

The item "/B*‘ is called a local sv;itch. It is local 
to the command line item "SC/^NER.BS" . ASM80 interprets 
a local B switch to specify the name of the output binary 

file. 

2*2.2 ASM80 Command Line Sv;itches 

A switch is a command line item which modifies 
the interpretation of the command line by ASM80. An 
example of a command line with sv/itches is: 

ASM80/L PGO TESTFOO/B (carriage return) 

In this example, "/L" is a switch on ASM80 commanding 
that a listing of the program to be produced on the de- 
fault listing device, "$LPT". The ”/L" is called a 
global switch since it modifies the command name. The 
switch "/B” is called local since it only modifies the 
processing of the filename is associated with in the 
cononand line, in this case **TESTFOO". 
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The following table lists the global switches 


that may be used by the v'rogranmer : 

Switch Function 

/K Change the default listing radix from hexidecinal 

to octal. If not present, all assembler output 

of program addresses and program data is in 
hexidecimal. This switch is equivalent to giving 
a "GRDX 8.” statement as the first statement 

of the program. 


/L By default no listing of the source program, 

and the corresponding binary output is given. 

- If a global /L is specified, a prcgram listing 
will be produced on the device "$LPT''4 


/N 


By default binary output is generated. If /'^ 
global is speciiiud, no binary output v.’ill be 
generated, formally AS.” SO v:ill generate the 
binary in a file created by appending ".38" 
to the name of the first program file in the 
Gommand line. 


The following table lists the local switches 
that may be used: 


Sv/itch 


Pun ct ion 


/E If no listing file is being output, this switch 

specifies the name of an error logging file. 

All error messages v;ill be dumped in this file. 

By default, if no listing is dumped and no /£ 
is given, all error messages will be dumped on 
the system console ($TTO) , If a listing file 
is specified, all error messages will go to 
the listing file, 

/L This switch specifies that the filename it speci- 

fies will be used as the listing file. This allov;s 
listings to files, other than the lineprinter. 


/B 


PMBlWAliPAmB 
e* FiXIB QQAIJTr 


By default a binary output is produced in a file 
created by appending ”.B3" to the filename of the 
first program file in the command line. If a 
/B is given, the filename it modifies will be ■ 
used fox dumping the binary output. 

These switches ate only those that are likely to 

be of use for the typical prograrrjner. Appendix 

gives the details for ether sv:itchcs used in testing 
and debugging the assembler* 
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2.3 Assembler Output 

■ASMPO oufont f.TilP lnt*n onr> of three classes* orrer 
mes s ages , 1 i s tings and bin ary . 

- t 

2.3.1 Error Messages 

When ASM80 detects an error in a program, 
an attempt is made to produce a reasonably specific 
error message to pin down the problem for the pro- 
grammer. Error messages are dimped on the listing 
file if one v/as requested or on the error file if a 
/E was given and no listing requested. Finally 
if no listing or error file v:as specified, the mes- 
sages are dumped on the system console. 

Error messages may be of three general classes ; 
Type Class 

E Error that v;ill almost certainly result in 

results other than those the prograr:ner 
desired. Examples v/ould be illegal register 
specification in an MOV instricution. 

W Warning of a possible error, but it is reason- 

able to expect the programmer might ’have in- 
tended this action. An example v;ouid be using 
this psuedo-op to back up the byte counter, 

D Diagnostic messace issued by the assembler to 

aid in debugging the assembler itself. Users 
should never get any diagnostic m.essages. 

In addition to these classes, some errors are 

of such a nature that it is impossible for the 

assembler to GOntinue processing the program. In . 

those cases the error is called Fatal, and will be 

Indicated as such in the error message. An example 

of a fatal error would be running out of symbol 


r.i.i ,-j -j '? K .H 




table space. 

The general format of error messages is: 
level nnrm error message text 
FATAt level nnnn text 

where level is one of the three level specifier 
characters (E,V7,D) , nnnn is the decimal error code 
assigned to this error message. If an error is 
fatal, it will be preceded with the v.'ork "FATAL". 

In the listing file, error messages precede 
the line of source they apply to. In the error’ file, 
the line of source on which the error v.’as detected 
is output also. 

In some cases, the error handler has a reason- 
able chance of identifying the actual source syirCjol 
that is primarily responsible for the error. In 
those cases , the symbol is included in the error 
message, enclosed in double quotation marks. In 
some cases, the enclosed symbol will be a 'blank — 
this generally means that the assembler v/as expect- 
ing a delimeter at this point other than a blank, 
i.e. , an operator in an expression. 

Appendix contains a full isting of all 

error messages that may be output by the assembler. 

■ » 

2.3.2 Listing Output 

If requested by the presence of a /L switch, 
the assembler v/ill produce a listing of the source 
program. This listing follows one of three formats 
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depending on the source line statement: 


• 

scatuiueitt 

OR 

iiiikkll 

statenant 

or 

nnnn iijjkkll 

statement 



In the first example, 

the statement required no 


processing by the assembler, i.e., a comment, thus 
the statement was simply reproduced in the listing. 

In the second case, the assembler had to evaluate 
an expression, but the statement did not generate 
ginary output, so the value of the expression is 
.listed, v;ith a precision of up to 4 bytes, j j , 

kk , and 11 each represent a single byte in hexideei- 
mal notation. Finally for the case of a statement 
that produces binary code, the byte address of the 
code is given by nnnn in hexidecimal, followed by 
the value up to four bytes in length. 

While the format will always follow this general 
organization, the precise field widths wi2».l be ad- 
justed to reflect the current output radix of the 
assembler. For example if octal output is selected 
the location field, nnnn requires 6 octal digits to 
represent any possible Intel address. 

2.3.3 Binary Output 

By default binary output is created for any , 
program assembled. This binary consists of a series 
of data blocks that specify the data to be loaded 
into the Intel and the address to load it. After 
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all the data blocks, an end block indicates the 
starting address of the program (see END psuedo-op) . 
Normally the binary will be produced on hexidecimal 
format, suitable for transmission over systera that 
can only process ASCII charaGters. For full details 
. of the Binary format see Appendix . 
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BUNK NOV BUiOi 


ABSTRACT 


The problems of conventional digital sampling are discussed 
in the context of patient monitoring and long term data collec- 
tion. Alternative methods of data compression are presented. 

The capability of three redundancy reduction algorithms to 
produce adequate representations of electrocardiographic data 
was examined. It was found that a zero order interpolator was 
the best redundancy reduction algorithm in terms of data 
reduction capability and overall quality of the reconstructed 
EK6. 

A hardware device was coi^structed which carried out zero 
order interpolation on a signal. Several alternative methods 
of implementing a zero order Interpolator were analyzed in 
terms of the demand they placed on a small general purpose 
computer. It was found that the eoE^uter could handle a moderate 
nundber of channels even with the worst implementation^, vdiile 
the best implementation produced negligable demands on the 
eonputer. 

Examination of the zero order interpolators reconstructed 
signal indicated that this representation was ad^Sjuate for 
analysis of rhythm. 
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CHAPTER 1 
INTRODUCTION 

Normal Heart Function (2) : 

The normal heart functions as two mechanically independent 
pumps. The right atrium (Figure 1) and ventricle pump blood 
through the lungs while the left atrium and ventricle pump 
blood into the systemic circulation. The mechanical action 
of the heart is coordinated by specialized tissue which 
forms the conduction system. This conduction system is rhyth- 
mically activated by a packemaker. 

The normal pacemaker is the sino-atrial (S-A) node (Figure 
1) . The S-A node initiates atrial excitation. As this 
excitation spreads, the atrial muscle is depolarized and 
• contracts , When the atrial excitation reaches the atroven- 
tricular (A-V) node (the only electrical connection between 
the atria and ventricles) it is delayed slightly before entering 
the Bundle of His. The Bundle of His bifurcates into the left 
and right branch, distributing the excitation to the Purkinje 
fibers in each ventricle. The Purkinje fibers interlace the 
ventricular muscle, distributing the excitation throughout the 
ventricles. The propagation velocity of the excitation through 
the Bundle of His (a / 4000 mm/second) is about an order of 
magnitude greater than the propagation velocity through the 
muscular tissue, thus all muscular activity is effectively 
synchrmized. 
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Figure 1 - Conduction System of the Heart. Drawing of the 
conduction system of the heart demonstrating the 
S-A node, A-V node, bundle of His, and the 
Purkinje system. After Corday and Irving ( 2) . 
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Muscular activity of the heart is reflected in the electro- 
cardiogram^ or EKG. The EKG is a recording of potential 
differences, vs time, between points on the subjects skin. 
Atrial depolarization (Figure 2) produces the P wave. While 
the cunplitude of the P wave will vaiy with the position of 
electrodes on the skin, the width is less than 0.10 seconds 
in the normal subject. Ventricular depolarization is reflected 
in the QBS complex. The width of the QRS complex is also less 
than O.iO seccmds. Pepolari ration of the ventricles produces 
the T wave. The Ta wave, generated by repolarization of the 
atria, is usually hidden by the QRS complex. The complete 
P-QRS-T complex occures over about 300 msec. 






Figuxre 2 
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- Excitation, and Contraction of the Heart During 
Normal Sinus Rhythm. The darkened area of the 
heart indicates the area the excitation spreads 
through to produce the darkened section of the 
EKG immediatly below. C the S-A node is always 
darkened) The first small deflection produced 
Ca-Atrial Exc.'.tation) is the P wave. This is 
followed by tlie large QRS deflection during 
ventricular eiccitation (d) , Finially the 
T wave is procuced by ventricular repolarisation 
during ventricular systole Ce). 

After Corday and Irving (2) . 
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iUbnormal Heart Function C ^ ) : 

Abnormalities of the heart may be classified as purely 
meGhanical as well as electrical. Electrical abnormalities 
often reflect structural damage; however, they may occure 
with no apparent mechanical damage. Abnormalities in the 
electrical activity of the heart may be classified as either 
disturbances of pacemaker function or conduction defects. 

Disturbances of packemaker function manifest themselves as 
disturbances of the S-A node or 'capture* of the pacemaker 
function by tissues other than the S-A node. An area which 
captures the pacemaker function (other than the S-A node) is 
called an ectopic focus. Since an ectopic focus initiates the 
excitation of the heart from an abnormal location, the se- 
quence of events in the heart cycle may be significantly altered 
cpreatly reducing the efficiency of the heart. Figure 3 
schematizes the sequence of events for one possible ventricu- 
lar ectopic beat. Such beats usually come earlier than eject- 
ed for a normal beat, and are called ventricular paremature 
beats (VFBs) . Figure 4a provides a stripchart recording of 
a VFB. The QRS complex for the VFB is premature and widened, 
due to the aberrant conduction of the ventricular excitation. 

Ohe specilized conduction system does not conduct the excita- 
tion as it would for a normal beat. 

The ectopic focus may be located in the atria also, in 
idileh ease the beat is called an atrial ectopic beat. The 
. excitation spreads abnormally through the atria resulting in 
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Figure 3 
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- Excitation and Gontraction of the Heart During 
a VPB, The darkened area of the heart indicates 
the area of the heaurt the excitation is spreading 
tiirou^ to produce the darkend section of the EK6 
below each heart. The first figure represents 
diastole of the previous (normal) beat. 

After Corday and Irving (2) • 
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Figure 4 
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Examples of a VPB,APB, and Nodal Beat. All these 
recordings came from the same patient. Note that 
his 'nonoal* P wave is abberantly conducted, 
resulting in a double peak. 

a) A VPB. The QltS complex is widdened, permature , 
and followed by a compensatory pause. 

b) A APB* Again the ectopic beat is premature. 

The QRS is normal, but proceeded by an abnormal 
P wave. 

el A nodal premature beat. The QBS is normal, 
but followed by an inverted P wave. 
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an abnormal P wave morphology. After the excitc^on is passed 
to the A-V node, the ventricular excitation sequence is normal, 
resulting in a normal QRS morphology. Atrial ectopic beats 
are also premature and are refered to as atrial premature 
beats (APBsl . Figure 4b provides an example of an APB. 

If the ectopic focus is in or near the A-V node, or the 
Bundle of His before the bifurcation, the resultant beat is 
termed a nodal premature beat. Atrial excitation progresses 
in a retrograde fashion, thus if a P wave is visable, it is 
inverted. The P wave may be buried in the QBS con^lex, which 
is essentially normal in shape. Figure 4c provides an example 
of a nodal premature beat. 

Conduction defects alter the normal passage of the excita- 
tion through the heart. The changes may take the form of a 
block, in which the conduction process is interrupted, or 
abnormal pathways. A block in the S-A node will result in 
failure of the heart to be paced by the S-A node; another 
focus must take over the pacing function. Typically the atrial 
tissue is most likely to take over the pacing function in the 
absence of S-A node initiation. A block at the A-V node pre- 
vents normal passage of the excitation from the atria' to the 
ventricles. Ventricuajj^ activity may now originate in the A-V 
node or in the ventricles themselves. However, independent of 
its origin, ventricular activity is now unrelated to atrial 
activity . 
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Clinical Significance of Ectopic Beats ( 2 ) : 

The clinical significance of the ryhthm defects mentioned 

in the previous section is not fully understood, VPBs are 

known to occure in otherwise normal subjects. Unless their 

severity is sufficient to cause patient distress (one can 

often feel the pause after a VPB) , nothing is usually done 

about them. The exact incidence of VPBs in normal subjects 

is not known. In the cardiac patient who is recovering from 

a heart attack, or other tramatic events involving the heart, 
t 

they are suspecjbd of being the first warning of impending 
heart failure through progressively more serioxis arrhythmias , 
such as Ventricular tachycardia and ventricular fibrilation. 
The frequency and variety (number of different foci) of VPBs 
is generally considered an indication of the irritability of 
the heart and influences decisions concerning the nature of 
the treatment given the patient. 






Patient Monitoring: 

Since VPBs,and other ectopic rhythms, are considered a 
measure of the irritability of a cardiac patients heart/ 
their detection by a monitor is of prime importance. A 
monitor for cariac patients should meet the following 
requirements : 

1) It should be able to detect ectopic beats in real- 
time and provide some record of their occurance, 

2) It should be operable even if the patient's 
'normal' rhythm, or beat morphology, would be 
considered abnormal, 

3) It should be insenstive to muscle noise, pacemaker 
and other artifacts. 

4) It may require training on a sample of the patients 
normal rhythm, but should not require training on 
ectopic beats. 

5) It should be as economical as is consistent with 
the above requirements . 

The first requirement, 1, is the most basic requirement of 
the monitor. The monitor must be able to monitor a patient 
continuously in real-time. Periodic monitoring of the patient 
offers little advantage over the current practice of intensive 
care units, where a nurse periodically monitors the patient. 

The ideal recoard of the occurance of an ectopic event would 
be a stripchart recording of the event and a few seconds of 
data around it. This record is in a form which the attending 
staff is familiar. Long term trend records would also be 
useful. 

Proper operation during a stable arrhythmia, requirement 
2, recognises the fact that any patient in a cardiac care unit^ 
in all probability^ has something wrong with his heart. 
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Thus his 'normal’ rhythm may in fact be abnormal for the 
general population. 

• M - 

Recognition of various artifacts that may occure in the 
EKG, requirement 3, is important. If possible the monitor 
should extract the EKG signal from the noise, but when this 
is not possible, it should be able to shut itself offi thus 
preventing false alarms. Excessive false alarms will result 
in inatteni^n to the monitor, a condition which must be 
avoided . 

The monitor must be able to detect previously unseen 
eVents, hence requirement 4 demands that only the patient^ 
normal rhythm be necessary for the establishment of detection 
criteria. 

Pinaiy, requirement 5 recognizes the rising costs of 
medical care in the United Sates. 

The author feels that the best approach to meeting these 
requirements, at the present time, is through the small 
general purpose computer. This is especially true at the 
research level. In Chapters II and III the author will 
examine one problem in the application of computers to 
the patient monitoring problem. 


Iiongterm Data Collection; 

As pointed out in section on the clinical significance 
of ectopic beats, the true incidence of ectopic beats in 
normal subjects is not known. One way to obtain this informa- 
tion is to record the EKG of a large population of subjects 
during normal daily activil^, then analyse the records. Say 
that one wishes to monitor a subject during normal daily act- 
for eight hours. If we assume that th^^lmportant compo- 
nents of the EKG Cin so far as the detection of ectopic beats 
is concerned) lie below SO Hz, and filter the waveform so that 
it is bandlimited to 50 Hz; then a sampling rate of 100 Hz is 
required aceordihg t& the Nyquist criterion C ^ * An eight 
bit sarnie provides adequate accuracy for this purpose, thus 

A 

2.3 X 10 bits are required to record eight hours of EKG. 
3here is simply no practical way to record this much data with 
out severely restricting the subjects mobility. Even if the 
data can be recorded, what does one do with it? Without some 
form of automated analysis, the systematic study of such large 
amounts of data is impossible at any reasonable cost. 

It Is apparent that a monitor meeting the requirements of 
the previous section would be suitable for processing such 
lon^erm records. If possible, it should operate faster than 
reatl-time to expedite the analysis of the records. 
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Conventional Digital Signal Representaion iS tO 1 : 

The conventio^nal approach to digital representation of 
signals is uniform sampling of the source signal and subsequent 
conversion of these samples into digital representation. The 
Nyquist criterion states that the sampling rate must be at 
least two times the period of the highest frequency component 
present in the source signal. If this condition is meet, the 
spurce signal may be recovered by lowpassing a uniformly spaced 
is^ulse train with amplitudes specified by the sample values. 

It should be noted that the Nyquist rate is determined by 
the highest frequency component present, not the highest 
one of interest. The usual technique is to reduce the band- 
width to the region of interest by lowpassing the signal, then 
Scunpllng the lowpassed signal at an appriate rate. In practice 
this does not produce tru^ly bandlimited signals, since all 
practical filters pass some components above their 'cutoff* 
frequency. Figure 5, after Macy(/o) , indicates the actual 
s^iqpling rates required for a specified accuracy , given a variety 
of filter roll-off rates and a non -bandlimited source signal. 

For this fiure the input is assumed to have a flat spectrum, 
and while this is not generally true in signals of interest, 
the figure is useful in that it indicates the magnitude of 
error one might e^ect. The possible error is quite large, 
and samples must be collected at higher rates than indi.cated 
by e^pling the Nyquist criterion to the cutoff frequency of 
the filter. 
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Figure 5 


(. 1 ) 


- Sampling Ratios Required For a Given Accuracy, 

As a Function of Roll-off Rata, The sampling 
ratio is (sampling frequency/filter cut-off freq- 
uency! , The errors are computed only for 
aliasing effects, errors due to nonoptimal 
reconstruction are not included. 

After Macy (10) . 
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While it might seem wise to have as sharp a cutoff as is 
possible, this may not be the case. In analog filters, 
caution must be used when 'simply' increasing the roll-off 
rate of a filter. It should be recalled that the phase 
characteristics of the filter are being altered. Poor phase 
characteristics will distort the data. In fact, the tin^ 
necessary for the transient portion of the filter fesponce to 
decay is inversly related to the 'sharpness* of the roll-off 
of a low pass filter. 

The selection of a sampling rate also depends on the nature 
of the data processing to be carried out on the digital repre- 
sentation of the signal. If the processing is to be carried 
out on a reconstructed signal, e.g, by digital lowpassing, 
then the signal may be reconstructed to an arbtrary degree of 
accuracy C limited only by the resolution offered by the number 
of bits used for the sample representation) . The errors 
presented in Figure 5 assume that optimal (lowpassing) recon- 
struction was used on the sample points. If optimal reconstruc- 
tion is not used the error will be even greater. 

If the application requires that the san^led data be presented 
for human observation, an additional constraint on the sampling 
rate must be applied. : .uman observers tend to interpolate data 
by connecting geometrica'.ly nearest neighbors. Thus a 59 Hz 
signal, sampled at 120 sauries per second, appears to the 
obsearver as two 1 Hz sinusoids in opposite phase. To guarantee 
that consecutive s angles are, in fact, geometric nearest neigh- 








bors requires a raidier high santpliug rate; soineUiiag on tlie 
order of 10 times the highest frequency component in the 
signal. 
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Alternative Approaches to the Representation of £KGs : 

The large amount of data that must be collected and stored 
for subsequent automated analysis of EKGs has resulted in 
interest in signal representation methods other than conven- 
tional sampling. Ideally these methods would be more efficient 
than sampling and easier to process. They must be practical 
to implement also. 

Young and Huggins (// } have represented the EKG with a sum 
of 12 empirically defined orthonormal exponentials. This repre- 
sentation is very compact, requiring only 12 values for each 
EKG complex (QRS and T wave only, they have not included the 
P wave in the representation) . Using these 12 exponentials 
they were able to represent a set of 18 EKG con^lexes with an 
error of only 5%. Iheir selection of exponentials was empirical 
and they make no claim that the set they have chosen can re- 
present any possible EKG complex to the same accuracy. 

The diagnostic capability of this representation has been 
studied (/s) . Young and Huggins found that in trials with about 
BO eos^lexes, they we able to separate the data into 5 disease 
classes and normals with about 70% accuracy. 

Cox et.ad. ( J) have proposed another technique for repre- 
senting the EKG. They call their method AZTEC, for Amplitude 
Z(»ie Time Epoch Coding, AZTEC reduces the number of bits 
required to describe an EKG by specif ing the waveform in terms 
of bounds and slopes rather than sample points. Bounds are 
specified as an^litudes and lengths (or durations). A bound 


provides the information that the input waveform was within 
a specified range of the bound amplitude for the bound length. 
During the rapid action of the QHS complex, many short bounds 
are produced. If a string of short bounds (< 8msec.) is 
encountered, and there has been no slope change in the input 
waveform, the string of bounds is replaced ’>y a slope segment. 
The slope segment specifies the waveform by a slope value 
and length. The input waveform may be reconstructed by inter- 
polating a straight line of specified duration euid slope from 
the amplitude specified by the preceeding bound. Cox et. al. 
C3) report that AZTEC is capable of providing a 10 to 1 reduc 
tion in the number of bits required to describe a typical EKG 


waveform. 
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CHAPTER XI 

REDUNDANCY REDUCTION 


Introduction : 

Figure 6 , after Kortman ( 8 ) r depicts various categories of 
data ccnnpression. All of these methods atten^t to reduce the 
bandwidth required to transmit a signal , that is the number 
of bits per unit time required to describe the signal to a 
given accuracy. These methods are parameter extraction, 
adaptive sampling, redundancy reduction, and encoding. 

Parameter extraction ( 8 Y reduces the number of bits required 
to describe a waveform by extracting the significant parameters 
of the signal and transmiting only the parameters . Thus if 
one were interested in the detection of VPBs one might con- 
struct hardwcure to extract the width of the QRS complex and 
the interval between co^lexes. 

Parameter extraction depends on a priori knowledge of the 
signal and a detailed knowledge of the purpose of the analysis 
to be preformed. In some cases the device to extract the 
parameters may be more complex than the analysis hardware. 

Adaptive sampling (8 ) provides bandwidth compression by 
adjusting the sampling rate of the signal to match the required 
bandwidth. Thus saiiples would be transmitted only when there 
is activity in the signal. 

Typically one is interested in g;^cu:anting that the digital 
i.'apresentaticm of the signal is accurate to a given error, 
rather than insist that it can be recovered perfectly. In 
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Figure 6 - Categories of Data Con^ression. This figure 

illustrates various irethods of data compression . 
After Kortman { 8 ) • 
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this case, redundancy reduction techniques iJ,6,g ) are 
useful. These techniques remove redundant information from 
the signal and transmit only that information required to 
represent the signal to a specified accuracy. 

Encoding (Q-) includes such techniques as delta modulation 
and incremental encoding. An incremental encoder transmits 
the difference between adjacent samples rather than the sample 
values. By adding the transmitted difference into a register 
(integrating the transmision values) one produces a representation 
of the signal at the receiver. Delta modulation is essentially 
a form of incremental encoding in which the transmitted values 
are limited to one bit. The receiver will add one unit to the 
register if a 1 is received and subtract a one if a zero is 
received . 
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Pr^di-Ctor?? and Interpolators : 

Redundancy reduction techniques fall into two general 
classes i 1 , Q\ , predictors and interpolators . 

Predictors make a predictions about the future behavior of 
the signal and transmit the parameters of the prediction. As 
the source signal varies in time, the actual value is compared 
with the predicted value. When the signal has deviated from 
the prediction by more than a specified amount, a new predic- 
tion is generated and transmitted. We will consider two 
types of predictors later in this chapter; the zero order 
predictor and the first order predictor. 

Interpolators generate a representation of a source signal 
based on the past behavior of the signal rather than predictions 
about its future behavior. Generally interpolators use some 
heuristic method for obtaining a good representation of the 
history of the signal. They may of course use an optimal 
technique, such as linear regression (if one is interested in 
representing the data with a line) ; but these methods’" are 
more dificult to implement. We will be concerned with the 
zero order interpolator, which is equivalent to the generation 
of bounds by AZTEC ( 3 ) . 
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Zero Order Predictor ; 

The zero order prdictor (ZOP) predicts that the waveform 
will remain constant in the future. Figure 7 depicts the 
operation of a ZOP. Consider the initial sample of the first 
sepient in Figure 7, at time t^. The ZOP transmits this sample 
value as its prediction and sets up an aperature of width 
2K centered on the transmitted value. This aperture defines 
the range the signal may occupy before a new prediction 
will be made, thus we always know the value of the source 
signal to with in tK. New predictions are transmitted at 
and t 2 since the signal left the previous aperture range 
on both sample points. The signal does not leave the aperture 
centered on the t 2 sample until time In this example, 

five values are transmitted while 12 sampljles were required. 

Under worst case conditions, note that the device would 
transmit every sample point, and thus become a conventional 
sampling system. 
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Figure 7 - Operation of a Zero order Predictor. The solid 

line represents the input signal, which is sampled 
uniformly in time at the points indicated by the 
circles. The squares indicate the transmitted 
samples. The reconstruced signal is represented 
by the line of alternating short and long dashes. 
Aperture boundaries are indicated by the dotted 
lines . 

i 

f 






38 


First Order Predictor; 

The first order predictor (FOP) iJ. ) approximates a signal 
with a line of specified starting amplitude, slope end length. 
Figure 8 depicts the operation of a FOP. Consider the sample 
at tg to be the initial sample of a new approximation to be 
gensrated by the FOP. The FOP waits for an additional san^Ie 
at t^, and then predicts that the signal can be represented 
by a line starting at the sample value of t^, with a slope 
specified by the difference between the san^le values of t^ 
and tj^. Apertures of 2K are centered on each predicted point, 
and the FOP checks to see if the signal sample value falls 
with in the aperture. As can be seen, the prediction fails 
for the sait^le value at t 2 » and a new prediction must be made. 

Note that the transmision of a prediction requires the 
transmision of two sample values, rather than one as required 
in the SOP. One might expect the FOP to be very sensitive 
to noise on the signal, since it bases the prediction on the 
derivative of the signal. In the werst case, the FOP'will 
limit out to a sample system^ since it must effectively 
transmit every sample point. 
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Figure 8 - Operation of a First Order Predictor. The solid 

line represent s the input signal which is unifomly 
sample in time at the points indicated by the circles. 
The squares indicate the values that must be 
transmitted. The reconstructed signal is indicated 
by the alternating short and long dashes. 
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Zero Order jLnterpolator : 

The zero order Interpolator (ZOI) { 6 , &) may be considered 

to be a heuristic for maximizing the amount of time the source 

signal remains within a floating aperture. Figure 9 depicts 

the operation of a ZOI. Consider time t^ to be the terminating 

point for a previ(^ segment of tlie approximation. We define 

3 , the maximum value of the source signal since the 

max 

start of a segment, and V_.^, the minmum value of the source 

mrn 

signal. V and V_. define the current boundaries of the 
max min 

floating aperture. Note, in Figure 9, that the aperture expands 
as samples for t 2 through t^^ are taken in. When the sample 
at tj ^2 taken in, the aperture expands beyond the maximum 
allowed width, K. An interpolated value is now transmitted. 

In the example of Figure 9 and in the results described 
later in this chapter, the interpolated value is the center 
'of the aperture at "" time timd at which the* error was * 

discovered. An alternative implementation of the ZOI, 
hereafter lefered to as a bounded ZOI , generates the interpo- 
lated value from the center of the aperature at the 

period previous to the detection of the out of bounds condition. 
In this case the errors of the ZOI are always bounded to i0.5K. 
The author chose the simpler, unbounded implementation since 
it does not require storage of the previous aperature center. 

In the worst case, that is the difference between san^le 
values exceeds K all the time, the unbounded system transmits 
half as many samples as there are samples. Each transrdtted 
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value is the average of two adjacent signal samples, thus the 
input signal is being digitally lowpassed. In the case of 
the bounded ZOI the device transmits every san^le point, 
and behaves as a conventional sampling system. 








o 
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Figure 9 - Operation of a Zero Order Interpolator. The solid 
line indicates the input signal, which is sampled 
uniformly in time at the points indicated by the 
circles. The dashed lines indicate the aperture 
boundaries. The reconstructed signal is indicated 
by the line of alternating short and long dashes. 
The squares indicate the transmitted values. Note 
that the reconstructed signal could not have been 
generated in real-time since the transmitted 
values occur e at the end of each segment. A device 
with a memery is required to reconstruct the 
signal. Since this memory must be bounded, the 
maximum lei:^th approimating segment is also bounded. 
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Results of E^qperimental Runs: 

In order to compare the ZOP,FOP, and ZOI on EKG data, 
the author has executed each algorithm on typical EKG data 
(.channel 5, PBBH-5 tape). Each algorithm was programed 
for the NOVA computer ( ) in the BASIC programming language 


C7 ) . Each algorithm was executed on exactly the same set 

of signal samples. The results are compared in Figure 10, 

a plot of experimental reduction ratio vs. per cent RMS error. 

The reduction ratio is defined as, 

number of output segments 
number of input samples 

The reduction ratio does not measure the actual reduction in 
bandwidth required to transmit the representation. 

The per cent RMS error is defined as follows; let the source 
signal valuer be represented by v®, vrtiere i runs from 1 to 
N, the total number, of samples. The reconstructed signal values 
for each sample point are given by V^, where i again runs from 
1 to N. Then the RMS Signal Value is given by. 



and the RMS Error is given by, 






A! 


Given the RMS Signal Value and the RMS Error, the Per Cent 


RMS Error is given by, 

RMS Error 
RMS Signal Value 
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Figure 10 


Plot of Reduction Ratio vs. % RMS Error. The 
reduction ratio and % RMS error as defined in 
the text. Each point in the polt was produced 
by executing a Z OP, FOP, and ZOI algorithm on 
the same set of input data. The value of 
K, the aperture parameter was varied from 10 to 
100. The reader should recall that for a 
aperture parameter of K, the ZOP and FOP set 
aperture widths of 2K, while the ZOI sets 
an aperture width of K. 







ft « ^ 
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Each data point in Figure 10 was generated by speciiing 
a value for the aperture parameter# K. The reader should 
recall that the aperture parameter specifies the maximum amount 
of allowed error before a new approximation is generated. Ten 
runs were conducted, over the same input data, with K ranging 
from 10 to 100 . The peak to peak signal range -was about 
600 units. 

Figures 11 through 21 illustrate the output of each 
algorithm for each value of K. The EKG data presented in these 
stripchart recordings is typical of the data used to produce 
Figure 10. For each figure, 11 - 21, the stripchart recordings 
presented are a) the original data, b) ZOI reconstruction, 
c) ZOP reconstruction, and d) POP reconstruction. 

Examination of the results of Figure 10, and the output 
presented in Figures 11 through 21 indicates that the ZOI is 
the best of the algorithms . The ZOI is least sensitive to 
the value of the aperture parameter and produces the most 
pleasing representation for all values oi K. 
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Eeconstructed Z0I,Z0P, and FOP Output for K = 10. 
This is the first of a series of figures for 
different values of the aperture parameter K, 

All the figures contain 4 stripchart recordings; 
a) the input signal, b) the ZOI output, C) the 
20P output, and d) the FOP output. The narrow 
spikes in the output of the FOP during the VPB 
are the result of a digital to anal§^ converter 
overranging, anf®^are not a result of the FOP 
algorithm. 





Figure 12 -- Reconstructed ZOI,ZOP, and FOP Output for K 







Figure 13 - Reconrstructeci Z0I,Z0P, and FOP Outpur for K = 3 
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Figure 15 - Reconstructed ZOI,ZOP, and FOP Output for K 














Figure 17 - Reconstructed ZOI, ZOP, and POP Output for K 
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Figure IB - Reconstructed Z0l,Z0P, and FOP Output for K 







Figure 19 - Heconstructed ZOI^ZOP, and FOP Output for K 
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CHAPTER III 

ALTERNATIVE IMPLEMENTATIONS OF THE 
ZERO ORDER INTERPOLATOR 

Hardware Implementation : 

A hardware system has been constructed which executes the 
zero order interpolation algorithm, A block diagram of the 
system appears in Figure 21. A brief discription of the 
hardware follows, for a detailed discussion the reader is 
refered to the Appendix. 

Refering to Figure 21, amplifier Q1 provides gain and off- 
set to condition the input signafor the following stages . This 
conditioned signal is sampled by the sample and hold amplifier 
Q2, which drives the maximum and minimum trackers. These track 
ers hold their current value until strobed by the TRACK com- 
raeind, at which time they aguire the new maximum or minimum 
voltage as required. The trackers may also be reset to the 
current sample voltage by the occurance of a RESET command. 
Amplifier Q3 computes the average of and , that is 

the current interpolated value for the aperture. Aii^lifier 
Q5 computes the current aperture width, an^lifier Q6 
compares this width against the maximum allowed width, Vj^. 

Thus a control level ERROR is provided when the aperture is 
too wide. Finally the Control Logic synchronizes the activity 
of the system with an external CLOCK aind provides* the required 
internal timing and gating functions , 


Figure 21 - Block Diagram of Hardware ZOI Implementation. 

This figure presents a block diagram of the 
hardware ZOI implemented by the author. This 
block diagram is functionally equivalent to 
the actual hardware; however, the reader is 
refered to the Appendix for a block diagram 
of the true system. 
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Figure 22 indicates the sequence of events in Uie system 

for a typical input siganl. Assume that the sample point s^ 

causes the system to terminate the previous segment and begin 

processing a new series of samples. Since BESET was high after 

S-, V and V . have taken on the value of S* . The leading 
0' max min 0 

edge of CLOCK initiates the next cycle starting with the acqui- 
sition of a new sample value, S^. Since the signal has increased 
in amplitude since the last sample, takes on the new value 

at the occurance of the TRACK command, continues to in- 

, max 

crease for samples and S^. When is taken in, the value 

of V . is decreased. This change produces an aperture width 
min 

which is too wide, and the ERROR level is raised. The occur- 
rence of an ERROR signal forces the sampling of V , thus 

»Vy 

generating the interpolated output value for the previous series 
of input samples. A RESET command is given after the next 
sample value , S^, is taken in. This resets the aperture 
width to zero. 

Figure 23 presents output of the system in the form of a 
stripchart recording of the direct output of Q4 . For the 
waveform presented, a triangle wave, the reconstructed signal 
and the output of Q4 are similar, however the reader is cautioned 
that they are not the same thing. A memory, to buffer segment 
<^nplitudes and lengths (which the hardware discribed does not 
include) ,is required to generate the reconstructed signal. 
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Figure 22 


(J ' 


Basic Hardware Timing Diagram, The reader is 
refered to the tejit for a discussion. 
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Figure 23 - Example of Direct Output of the Hardware, 

A triange wave of 0,5 Hz, is the input 
signal Ca) ^ the ouput is presented in (b) . 
Tiie reader is reminded that this is the 
direct output not the reconstructed output. 




' ' r-rp r 




68 


LJ Processor Overhead for Alternative Implementations 

of the Zero Order Interpolator; 

In the context of processing the output of a ZOI on 
a small computer , the percentage of central processor 
unit CCPU) time that must be devoted to execution of the 
ZOI algorithm is iit^ortant. One would like to make as much 
time as is possible availiable to analysis of the data, rather 
than just compressing input data. 

We will consider the percentage of CPU time required for 
three possible implementations of the ZOI algorithm. The 

V 

following three implementations will be considered; 

I. Direct software controlled sampling of each signal 
channel with a software implemented ZOI processor. 

II, Specialized hardware controlled sampling of each 
i signal channel with menySory access through a direct 

memory port and a software implemented ZOI processor. 

Ill, Specialized hardware implementation of a ZOI 

processor with software controlled transfers of 
segment amplitudes. 

In addition, we will consider the cost of conventional sampl- 
ing as a reference point. 

Table 1 presents the execution time of various software 
processes required to implement the above systems. The times 
cire based on assembly language programs written by the author 
and the timing information provided in 14 -) for the NOVA 
1200. 

In addition to the timing information of Table 1, some 
assumptions about the nature of the source signal must be 
I made inorder to estimate the preformance of the software and 






Table 1 - Software Event Timing 


Event Time to: process event (/isec) 


a) Interrupt from any external 

device 40 

b) Input sample changes Vmax 37 

c) Input sample changes V^in 40 

d) Input sample changes Vmax and 

puts aperture out of bounds 49 

e) Input sample changes Vmin and 

puts aperture out of boiinds 51 

f) Buffering output points of ZOI 30 

gl Identification of channel number 

in hardware zor 11 

h) One data transfer by direct memory port 1,2 


All times are based on assembly language programs 
written by the author. Times are execution time 
on a NOVA 1200 i4) , 
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and hardware. These assumptions are; 

I) A randomly selected sample point is equally likely 

to change as V„ . . 

max min 

II) Every input sample will change either V or V_^ 


max 


min 


Assumptoin I is straightforward. It should be noted that 
assumption I implies that a randomly selected ZOl output 
segment is equally likely to teminate on a change in as 


min 


. Assumption XI is, in fact, not true; but is made to 


simplify the estimation of expected processing time per input 

point. It can be justified on the grounds that the software 

requires about as much time to process a point that does not 

change V or V . as it does to process a point that does 
max min 

change one of them. 

A general purpose computer with analog to digital conversion 
equipment and a analog channel multiplexer forms the hardware 
of implementation I, The software costs can be broken down 
into three arl|^s, (1) interrupt overhead to divert the processor 
from any other task it might be preforming, (2) the expected 
processing time in the ZOI software for each sample point, and 
(3) the expected buffering time to save any output from the 
ZOI software. The time required to respond to an interrupts 
fixed at 40 //sec. The actual time required in the ZOI soft- 
ware is a random variable, however using the assumptions > 
presented earlier, and specifing a reduction ratio, R (see 
the previous chapter) , we may compute the expected value of 
the processor time required per input sample to be. 


71 


40 + C38.5 + 41.5R1N microseconds, 
where N is Uie number of active channels that must be compress- 
ed. This expected value may be converted in to peiQcentage 
of CPU time by dividing by the period between input samples. 

System II uses some what more complex hardware in an 
efort to avoid the relatively high interrupt overhead (for 
small numbers of channels) by placing the sampled data in 
memory through a direct memory port. 'Hiis avoid the interrupt 
overhead for each point sampled, since an interrupt would not 
be generated until the core buffer was full. One d^advantage 
of this method is that it does require buffering the original 
data and thus req^iires core that Systems I euid III would not 
need. Refering to the timing information of Table 1 and the 
assumptions of this section, the expected value of the 
processor time required to handle each sample point is, 

C 39,7 + 46.5r)N microseconds, where 
where R is the reduction ratio, and. N is the number of active 
channels. As can be seen, this method offers little improve- 
ment over the previous system. 

A specialized hardware in^>lementation , as. described inthe 
previous section, comprises system III. In this case the 
processor time per output of the ZOI hardware is fixed at 
81 microseconds (interrupt overhead, buffering, and channel 
Identification) per outpu'^oint. Since the fraction of s^le 
^points expected to result in output of a segement is given 
by R, the expected value of processor time per input sample 
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O point is, 

81PN microseconds , 

where R is the redundancy reduction ratio, and N is the number 
of active cheinnels. 

The alternative systems are compared in Figure 25 
for R = 0,1 and a basic sample interval of 2000 microseconds 
C500 Hz.l. The plot presents percentage CPU time vs, the 
number of active channels. In addition to the three systems 
discussed in detail, the overhead for conventicj^l sampling is 
presented . 

Both implementations involving software ZOI processing 
are comparable to conventional sampling, while the hardware 
implementation of the system is superiojrfr. 

( I' ' - 
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Figure 24 


t 


Per Centage of CPU time vs. Number of active 
Channels. The plot is for an assumed reduction 
ratio of 0.1 and a baUgc sampling rate of 
500 samples per second. The per cent axis may 
be scaled for other sampling rates. For instance 
for a sampling rate of 250 samples per second, 
the value of 2C% becomes 10%. 
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CHAPTER IV 
CONCLUSIONS 

General : 

The object of this work was to determine an efficient and 

practical representation of the EKG for digital processing. 

The results presented in Chapters II and III indicate that 

the zero order interpolator provides an efficient representa- 
. ^ 

tion of the EKG and xs practic^ to xmplement,both in soft- 
ware and hardware. 

Examination of the records presented in Figures 11 through 
21 indicates that for moderate values of K (around 30) the 
ZOI provides an excellent representation of the original signal. 
Certainly no cardiologist would have any difficulty making 
a rhythm analysis from the results presented. Redundancy 
reduction ratios of 0,1 to 0.05 fr(g moderate values of K indi- 
cates that the representation is indeed more efficient than 
conventional sampling. ^ 

The results cf Chapter III, in particular the data siimmarized 
in Figxire 25, indicates that the implementation of the ZOI 
algorithm in software places a burden on a processor roughly 
equivalent to the demands of conventional sampling. If the 
algorithm is executed in special hardware the processor demands 
are negli gable for any reasonable number of channels. 


J 
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APPENDIX 

HARDWARE DISCUSSION 
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General Information; 

The hardware design for the zero order interpolation 
algorithm uses both aoialog and digital integrated circuits. 

The analog circuitry is constructed on two Vector boards 
(Vector Electronics Co., Inc.), while the digital logic is 
constructed on CBG-5 logic boards. The CBG-5 boards have 
space for three 14 or 16 pin dual in-line packages. The 
edge connection is 44 readout, and all package pins other 
than power and ground are available at the edge connector. 

The cards were assembled in a card rack, and all intercon- 
nections are via AMP Taper Pins (AMP Inc.), The analog 
functions were implemented with operational amplifiers; 
the LM301A differential amplifier, the LM310 follower, eind 
the LM311 comparator (National Semiconductor Corp.). These 
amplifiers were chS^en for their high preformance/cost 
ratio. The control logic was implemented with series 74 
transistor- transistor logic (Texas Instruments, Inc.), This 
logic was chosen because of its availability and the medium 
scale integration functions, such as counters and decoders, 
included in the series . 

Figure A1 provides a functional schematic of the system. 

The hardware divides into two major sections, the analog 
circuits and the control logic. The analog circuits provide 
functions such as sample and holding, maximum tracking, min- 
imum tracking , and voltage comparison. The control logic 
accepts inputs from an external CLOCK and the analog circuitry 
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and generates coinmands to the analog circuitry. 
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System Control Logic: 

The control logic is presented in Figure A2 and an 
associated timing diagram is presenred in Figure A3. The 
control logic generates a series of time states (TSO ,TS1 , . . , , 

TS9) which are gated with other logic levels to provide 
coiifmemids to the analog circuitry. 

Assume the system is waiting for the next CLOCK signal 
to start a new cycle. Refering to Figure A2 , consider element 
El t a D type flip-flop) to be in the low state CQ=0) . Its 
Q output thus asserts the reset function of element E4 (a 4 
bit counter) . Since counter E4 isreset, all its outputs are 
low and element E5 ( a decoder) decodes this as a zero. The 
ten decoded outputs of E5 correspond to the time states of the 
system, the the system is now in TSO (no output is shown 
on E5 for TSO since this logic level is not actually used by 
the system) . When a zero to one transition occures in the 
CLOCK level, El changes to high state. This enables the 
counter, E4, 

Elements E2 (an open collector nand gate) and E3 ( a Schmitt 
trigger) form a gated multivibrator which clocks the system. 
Ignoring E2 for the moment, we will explain the operation of 
E3 and its associated passive componenets. When power is first 
turned on, C is discharged, thus the input of E3 is low, forcing 
the output high 4. 5 v) . When the output goes high, C begins 
to charge up, to 4.5v, through resistor R2. Note that diode 
d is backbaised so that Rl doe^ot affect the charge up cycle 
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of C. Wh.en C reaches the upper threshold of the schmitt 
trigger, the output of the trigger will go to logic 0 (/>^0.4v). 
Now C begins to discharge through both Rl and R2 until it 
reaches the lower threshold of the schmitt trigger. When 
this happens, the output will go high , and the cycle will 
then repeat. Element E2 gates the muti vibrator by grounding 
tht,e input of E3 when both of its inputs are high. Since 
input E2-1 is from the output of E3, the multivibrator can 
only be shutdown when E3 is in the high state. This forces 
an integral number of puloes of constant width to be produced, 
with out reguard to the relative time the gating input, E2-2, 
is raised. The values of R1,R2 and c where chosen to provide 
a clock period of about 20 microseconds, thus each timestate 
(except TSO) is 20 microseconds in duration. 

The clock pulses produced by E3 are counted by E4, and 
E5 decodes the appropriate time state. When TS9 is decoded, 

El is cleared through the asynchronous clear, pin 1. This 
shuts off the multivibrator and clears the counter causing 
TSO to be decoded. The control control logic will now wait 
for another external CLOCK to initiate the next cycle. 

Time states Tl, TS3 and TS5 are gated to produce commands 
to the analcc! circuits. Refering to Figure A3, note that 
ERROR is high at the start of cycle 0. Thus the system error, 
or aperture width, went out of bound during the previous 
cycle. During TS5 of the previous cycle, E9 was clocked 
C Figure A2) into the high state since ERROR was high. 
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Cycle 0 starts with the CLOCK transition, and TSl is 

entered. TEl always generates an INPUT SAMPLE command 

since the input sig*^! is sampled every clock cycle. E9 

is in the high state from the previous cycle, enabling 

E12, This results in an OUTPUT SAMPLE command d, :.*.ng TSl. 

OUTPUT SAMPLE forces the output sample and hold to acquire 

the current value of V as the output of the system, V , 

avg out 

^out represents the interpolated value of the 

input signal since the last OUTPUT SAMPLE command. The inter- 
val between OUTPUT SAMPLE commands is the length of the seg- 
ment . 

TS3 is entered next. Since E9 is in the high state, 

ElO generates a RESET command. RESET forces the Maximum and 

Minimum trackers to acquire the current value of the input 

sample, S^ in this case, V and V„. are at the same value 

^ 0 max min 

so ERROR retuams to the low state. TS5 will now clock E9 into 
the low state. 

It 

During cycle 1, The system acquires a new input sample, S^, 
during TSl, No OUTPUT SAMPLE command is generated since E9 
is low; however, TS3 generates a TRACK command through Ell, 

For this cycle, V„ increases in responce to the TRACK com- 
mand. ERROR remains low, and L3 is clocked low by TS5, This 
sequence of events repeats for cycles 2 and 3. 

^min changed during cycle 4, and the aperture width 
increases beyond the limit, Vj^, forcing ERROR high. TS5 will 
clock E9 high, ths(u forcing ah OUTPUT SAMPLE command in cycle 5. 
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Analog Design: 

Signal conditioning for the system is provided by amplifier 
Q1 fo Figure Al, Figure A4 is the schematic diagram of this 
stage. The stage is of conventional negative feedback design, 
providing a gain range of 0.1 to 2. An adjustable offset vol- 
tage is applied to the positive input of the amplifier through 
the resistive divider of R5 and R6 , providing a DC offset of 
±5 volts. Cl provides feedforward compensation for the ampli- 
fier, as recomended by the manufacturer. This increases the 
slew speed of the amplifier to 10 V//xsec., resulting in im- 
proved bandwidth and step responce. A wideband amplifier is 
not required in this stage and a less expensive amplifier 
could be subsituted. 

The analog sample and hold function was required to sample 

the input voltage Camplifier Q2 of Figure Al) and to hold the 

output of the system (amplifier Q4) . In addition the Maxi^aum 

and Minimum Trackers use this function. In all cases the 

sample and hold function was implemen^,ed with the circuit if 

Figure A5. On command the storage capacitor. Cl, is charged 

to the value of the input voltage, switch, 

Ql. The LM310 voltage follower allows readout of hte stored 

voltage. Pertinent specifications of the LM310 follower and 

LS4391 FET switch are: 

LM310 Voltage follower 

Input Bias Current 
Input Resistance 
Offset Voltage 


10 . ^ (maximum) 

10 ohms (minimum) 
10 mV (maximum) 






s>^ 



88 


LS4391 (2N4391) FET switch- 

R on 30 ohms (maximujn) 

Off state drain leakage 0.1 nA (maximum) 

Pinchoff voltage 10 V (maxi mum) 

Gate-source capacitance 3,5 pf (maximum) 

TSie major considerations in a sample and hold circuit 
are the acquisition time and droop. Acquisition time is the 
time required for the sample and hold to come to within a 
specified error range of the input signal, when switched into 
track C or sample) mode. Droop is the decay rate of the out- 
put voltage while in hold mode. 

The FET and capacitor, of Figure A5, form an RC circuit, 

■6 

In the on st^e the worst case resistance of the FET is 30 
ohms, thus the RC network has a time constant of 0.3 micro- 
seconds, with Cl = 0,01 microfarads. Tliis circuit will require 
about 4,6 time constants to charge to within 1% of l5^e final 
value, thus the capacitor acquires the input signal to within 
1% accuracy in about 13.8 microseconds. The LM310, due to its 
high slew speed and excellent lar^ signal step res^nce adds 
little to the overall acquisition time of the sample and hold. 

Droop is caused by charge leaking off the storage capacitor. 
This charge may leak out through the FET switch and the input 
of the LM310 (provided good capacitors are used the internal 
leakage of the capacitor is not important) . The major source 
of leakage current is the bias current of the LM310 . Given 
a worst case bias current of 10 nanoamperes , the worst case 
voltage droop is given by. 
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Thus in 100 msec., a reasonable maximum holding time, an error 

of J7(r mV nay accumulate. For a peak to peak range of 5 volts 

3 . 

this would correspond to error. 

The FET is switched by transistors Q2 (Figure A5) and Q3 
through diode Dl. When S^iMPLE is high, Q3 is on and its col- 
lector is at ground. This pulls the base of Q2 down through 
R5(, biasing Q2 on. When Q2 is on, its collector is at +15 
volts, thus backbiasin^ diode Dl. Sine the drain and source of 
the FET are connected through Rl, they are now at approximatly 
the same voltage and the FET is in its on state. 

\flhen SAMPLEis low, Q2 and Q3 are off. The collector of Q2 
is pulled to -*15 volts through R3. Diode Dl is foward biased 
for any V. above -15 volts, thus the FET ivuiction is 
reverse biased. Pinch off voltage for the LS4391 is 10 volts, 

’i. 

meucimum, bo that input voltages above -5 volts will keep the 
FET in pinchoff. 

Considering just the PET, the turn on time is given by the 
time constant of Rl and the drain-source capacitance (3.5 pf) 
of the FET, or 0,35 microseconds. The turn off time is deter- 
mined by the timeconstant of R2+R3 and the drain -source capaci- 
tance, or 21 nanoseconds. Both of these times are small compared 
to hte acquisition time and can be ignored in this application. 

The maximum tracking function was implemented with the cir- 
cuit of Figure A6 . This maximum tracker is not a continuous 
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tracker, that is it can track the changes in the input only 
when gated by the TRACK command. In addition the tracker 
may be reset at any time to the current input value through 
the RESET command . 

A major part of the maximum tracker is a sample and hold 
circuit which was discussed earlier. The circuit for the 
maximum tracker is presented in Figure A6 . The LM311 compara- 
tor, compares the input value to the tracker, with the 

current output, If the input is greater than the output, 

the output transistor of the LM311 will be off, thus the 
base of Q5 will be pulled up by R9. Q4 and Q5 form an *and* 
gate. If is greater than (implies Q5 is on) and the 

TRACK command is high (thus switching Q4 on) then the base of 
Q2 will be pulled to ground through R5 , Q4 qnd Q5. This will 
switch Q2 on and thus switch the FET on^ 

Input offset adjustments (rll,R12,R13) have been used on 

both the IiM3I0 and LM311. Note that at least one of these 

offset adjustments is required. Assume that the LM310 is ideal 

and has a zero offset voltage, but the LM311 selected for the 

circuit has a wor^t case offset voltage of +10 mV. If and 

Vout were eactly the same value, the LM311 would say that 

was less than V. all the tme. In this situati0on the system 

in 

would always track the input signal at the TRACK coimnand. 

To achieve the minimum tracking function the unity gain 
inverter of Figure A7 was placed in front of another maximum 
tracking circuit. The output of this minimum tracker is the 
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seojn reasonable to build a minimum tracker by simply reversing 
the input to the LM311 in the maximum tracking circuit. This 
method will not work, however, due to the positive input 
current of the LM310 follower. Since the LM310 has a positive 
input bias current, the output of the maximum tracker will 
slowly decay with time while in the hold mode. Consider 
the situation of a constant input voltage (input to the ZOI 
system) with a small amount of noise on it. In this situation 
the maximum tracker would take on the value of some local 
maximum and begin to decay. When it had decayed to much, it 
would be reset by the LM311, since it must eventually decay 
to a level below the signal level. If we simply switched 
the inputs to he LM311, this circuit would not work as a 
minimum tracker since it would be decaying away from the 
value of the signal. Since it will always be below the signjlal 
once it start to decay, it will never be reset to a value 


close to the true signal value. 

Amplifier Q3 of Figure Al computes the average value of 

V and V . , or the center of the aperature. Since -V . 
max min nun 

avcdlable from the minimum tracker, one wishes to compute 

0.5 (V -C-V . )), This is accomplished with the circuit of 
mxri 

Figure AS, a differentail configuration. If R3 is adjusted 
so that the ratios Rl/(R2+R3) £ind R4/R5 are eq^yl, the -ther 
output of the stage is given by. 
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The ratio R4/R5 = 0.575, which is equivalent to a gain 
change in the signal. 

Feedforward compensation is used in the averaging stage 
to improve the step responce of the LM301A. Thi si done 
so that the output of the stage may be used in the TS follow- 
ing the change at the inputs , 

The error comparator ,Q5 of Figure Al, provides the ERROR 
signal if the aperture width exceed the threshold Vj^, Figure 
A9 provides a detailed schematic of this circuit. The pertinent 
LM311 input specifications are; 


Input offset 

voltage 

10 

MV 

(maximum) 

Input offset 

Cureent 

70 

nA 

(maximum) 

Input bias current 

300 

nA 

(maximum) 


The low bias currents of the LM311 allow the resistor 

network of R4 and R5 to be used to compute the difference 

V -V . . Assuming R4 and R5 are driven by ideal voltage 
max min 

sources, the voltage at the junction of R4 and R5 is 
0.5 CV -V . ) . The LM311 compares this with the reference 
voltage, Vj^. Rl, R2, and R3 provide a reference votage over 
the range 0 to 0.5 volts. Cl bypasses noise at the minius 
input of the LM311, Refering to Figure A9, note there is 
a 1.2 K ohm inbalance in the input impedance to the LM311, 
however, even under wort case conditions, the low input and 
bias currents of the LM311 hold the maximum error to 0.5 mV. 
This error can safely be ignored. 

The output of the LM311 (pin 7) is an uncommited collector. 
This collector is tied to the junction of R7 and R8 , R7 and R8 















' '—f '^I~r tfinUlT^tlWii*' 1 I ■ ■• j fiTI >1- • r i"r " 






bff A Vi'.' l-tfifeal'E.'w a ’r. -^1 


IN 

o 




98 


Hardware Syatein Output: 

Direct output of the hardware system is presented in 
Figure AlO, Figure AlO-a shows the output for a triangle 
wave input of low frequency. It is stressed that this is the 
direct output of the system, not the reconstructed output. 

A device with memory is required to generate the reconstructed 
suganl representation, and the hardware system described in 
this appendix has no such memory. 

The steps in the direct output occure with the OUTPUT SAMPLE 
command. Thus a given segment in these figures provides amplitude 
information for the previous segment length. 

jr 

Figure AlO-b provides an example of the ^stem output 
with a triangle wave of increasing amplitude. This gives the 
reader an idea of the nature of the representation for various 
amplitude signals, relative to the aperture width. 





Fijure ^JO o-f h)r^ 0».ipoJ- 
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ABSTRACT 


A system to analyze magnetic tapes of electrocardio- 
grams i£CG*sJ recorded over 2^ hour intervals has teen 
designed. Tapes are examined on a oeat-to-teat Oasis for 
potentially significant arrhythmias at a rate substantially 
faster than required to acquire the data. 

Using a L>ata General hova computer and a special 
design interface system, the magnetic tape recordings are 
sampled 200 times per second ^real-time) . This amount of 
data is reauced by a factor of io by creating a horizontal 
line segment approximation. This is followed by a block 
idiieh detects the QR3 complexes, another which measures 
the width »T\d magnitude of the QRS complex plus the baseline 
and QQ interval, and another block which uses these measure- 
ments to decide if the complex represents a FVC or a 
normal heartbeat. Output consists of coded representations 
of the width, magnitude, and QQ inlerval, plus a flag 
marking detected PVC's. 

This system has been tested on several tape recordings 
of clinical data, containing over 75 premature ventricular 
contractions (PVC*s>, Every PVC was flagged. The rate 
of false positives is under with the expectation of 
being greatly reduced. The program runs at 8 times real- 
time, calculations indicate that it presently should 
•be able to run times real-time on the Nova. It is 
ei^eeted that the same program should z^in 60 times real- 
time on a Nova 800. Reducing the sampling rate to 120 
samples per second is probably feasible, which would yield 
a corresponding increase in speed. 
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BACKGROUND, AND THE PROBLEM TO BE SOLVED 


The electrocardiogram (ECG) la the potential created 
on the surface of the body by the electrical activity of 
the heart. Because many problems with the heart are in- 
volved with a defect in its electrical generation or con- 
duction apparatus, the ECG is a very useful clinical tool. 

"No cardiac study is complete without an electrocardiogram. 

No electrocardiographic analysis is more detailed than a 

20 

computer aided one." 

The Normal Heart 

The heart is composed primarily of several layers of 
muscle. When they contract in orderly fashion, blood 
is circulated to the entire body. The rhythmic contraction 
is controlled by a specialized electric generation and con- 
duction system in the heart, as seen in i.igure 1. An impulse 
is started each beat in a natural pacenaher in the right 
atrium called the sinoatrial node. The impulse travels 
through the atrium, depolarizing it and causing It to con- 
tract, and stimulates the atrioventricular node. It ontin- 
ues along a comparatively high speed system calledL the Purkinje 
fibers, spreading through both ventricles and causing them 
to depolarize and contract. After a short period of time, 
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Figure 1 


Physiology of the heart | including pacemaker 
and conduction system 
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the heart repolarii&tis and la ready Tor aiiotheir cycle. 

The movement of the electrical Impulee through the heart 
can be reasonably modelled as an electric dipole^ It 
obviously has a three-dimensional pattern of movement. This 
equivalent dipole is observed on the body surface by two or 
more electrodes, so that the placement of the electrodes 
determines which planar projection of the dipole is observed. 

The observed waveform is dependent on skin resistance, 
amount of body fat, heart placement, and other factors in 
addition to electrode placement, so that it becomes impor- 
tant to locate the electrodes skillfully to obtain a useful 
ECG. Figure Z depicts a stereotyped £CG. The P-wave 
corresponds tv the atrial depolarisation, the QRS complex 
to the ventricular depolarisation, and the T-wave to ventri- 
cular repolarisation. 

Arrythmias 

Problems with the hearths electrical system may be due 
to the impulse origin or to abnormalities of the pathway, 
and may be either regular or intermittent. Many classifi- 
oatlons nave been described, and one particular one la the 
premature ventricular contraction (PVC), 

In a PVC, the focus, or impulse origin, is not In the 
sinoatrial node, but in a ventricle. It Is earlier In a 
cycle than the normal Impulse would be expected, and, therefore, 
overrrldes it, it is generally unable to make use of the 
Furklnje fiber conduction system, and consequently travels 
slower and creates a wider QRS complex. Figure 3 shows three 
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Figure 2 

Normal ECG waveform 

P-wavei atrial depolarization 

QRS Gomplexi ventricular depolarization 

T-«avet ventricular repolarlzatiom 


Figure 3i Three PVC's* illustrated In the context 
ef their normal heartbeats. Note the prematurity 
In all of them and the varying degree of 
transformation of the QRS morphology. 
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different FVC*s in the context of beats which are noxinal 
for that patient. If all FVC*s have the same ventricular 
focus, they vdll tend to have the same shape. If the FVC's 
are multi- focal, they may h£ive widely varying morphology, 
and this Is syaptonatic of different problems than unifocal 
FVC*s. The FVC is always early, usually has a compensatory 
pause following it, is usually wide, is often taller than 
normals and inverted in polarity, and may have a bisarre 
shape. 

Clinical Significance— Research 
Isolated FVC*s may occur in otherwise normal individ- 
uals, and probably have littly significance as long as they 
are infrequent. As they become more frequent, and if they 
are multi-focal, the danger is increased. In partictilar, 
if early FVC*s fall on T -waves (the electrically vulnerable 
period) they are likely to initiate fibrillation, a dis- 
organised beating of the heart which has very little effec- 
tive pumping action. Fibrillation is reversible if caught 
within about three minutes, else death usually results. 

Until recently, the patient depended on very fast recognition 
and response to his fibrillation for his life. Since 1965# 
in some coronary care units, parients with FVC*s are treated 
with anti-aXThythmie drugs such as lidocalne to suppress 
FVC*s and prevent fibrillation altogether, Frior eaq)erxence 
had shown that about 159^ of heart attack patients in the 
horpltal developed fibrillation. Since then, incidence of 
vantrieulmr fibrillation has been held to less than 
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This technique has reduced the mortality rate of those in 
coronary care units from 3o;2 to 20?S, or about 60,000 lives 
per year. This is certainly sufficient to establish a FVC 
monitor as an important clinical resource. 

In order to do scientific research on such topics as 
testing anti -arrhythmic drugs, l::;;rge numbers of patients 
must be monitored for long intervals. Another research 
area of critical importance is to attempt to find predis- 
posing symptoms in those who ** drop dead** from sudden heart 
attacks. Towards this end an entire community has their 
EGG *8 recorded for a number of years. There is currently 
a ioint U.S. -Soviet program underway to study thousands of 
people.^ Here again Is a huge mass of data to be analysed. 
Because PVC*s arc generally not regular with respect to 
time, a large sample of data must always be ar^alyzed. There- 
fore, a method for making an accurate, high-speed detector 
of FVC *8 is urgently needed. Human observors are too 
scarce, expensive, and inaccurate. 
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i^HfivIOuS WuHK 

There has been a great deal of Interest in applying 
computers to the analysis and monitoring of BCG's, siroe 
about i960. Every researcher seems to have his own special 
idea, consequently there has been a proliferation of hard- 
ware used, detection methods, end decision algorithyis. 

Many have been quite successful for either general analysis 
or a speciaillaed job they have been designed for. One 
excellent source that covers the entire field f^om many 
perspectives is a 1970 book edited by Cesar K Caceres and 
Leonard S. Dreifus.^ Seventy-one authors have contributed 
articles ont instrumentation, techniques In computer pro- 
gramming, utilization and potential, and research for clin- 
ical needs. The book covers both the state-of-the-art and 
a systematic discussion of needs for the future. There is 
one paper on monitoring, none on high speed analysis. A 
1972 paper by Jerome R. Cox and colleagues has a fine review 

an^ iiscussion of who is doing what in computer analysis of 

6 

BC0*s, EEC *8, and blood pressure waves. His paper Is also 
notable for the 300 references included on the field. Many 
papers on rhythm monitoring are detailed, but, again, none 
on hlgh-spoed evaluation. 

High-speed ar^alysis, ev m specialized to a single feattire 
such as yVC*8, Is not described in the literature. One technique 
Is to display magnetic tape recordings at 60 times real-time 
on an oscilloscope that triggers such that normal waves are 
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BUperimposed. Abnormal waves stand out vlsuaJ.iy to tne 
human operator, and the system may be halted at each FVC 
to provide a strip chart recording. Heart rate may be mon- 
itored as an added feature,^ One computer analysis runs 

12 

as fast as 15 times real-time on a PDF -7, Simple analog 

devices such as high and low rate alarms may also be used. 
This author is unaware of other working systems in high- 
speed evaluation. 

The data reduction scheme used in this paper, which 
is a critical part of the program, is based on a paper by 
J.B.WaXters, Jr,^^and is similar to an algorithm used for 

7 

a similar purpose by Jerome R.Cox, Jr.' 
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THE SOLUTION 


Kequiretnents 

Firsts the analysis must he accurate in its decisions* 

The numher of false negatives, PYC*s classified as normals, must 
he very close to zero, and the numher of false positives, normals 
' and artifacts classified as ahnormals« must he reasonably small. 

The acceptable number of false positives depends on the respective 
costs of false negatives and positives, if the analysis Isto he 
used clinically or for research, and on the condition of the re- 
cording, On a patient with a heart rate of 100 heats per minute 
and 1 PVC per minute, a O.l^ false positive rate will have 1 
false positive for every 10 true positives and a total of 86 
false positives In a hour recordizig. This would probably he 
a good level. The numher of PVC*s is much smaller than the 
numher of normals, and it is they that hold the information that 
may indicate disease in the patient. Also, if the patient has 
more than one type of PVC , it is very important not to catch some 
types and miss others* This is because runs and p^-ttems of PVC*s, 
and the realization that they are multi-focal, can be as critical 
as finding individual PVC*s. Doctors may also he particularly 
prone to doubting the validity of a program that misses an arrhy- 
thmia that they spot visually, and without the physicians coop- 
eration, no automation will be used. **Hours of accurate 
detection of normal QRS complexes may he less important than 
the identification of a few anomalous heats that possibly 
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4 

warn of more severe arrhythiulae to couie.** Thr so w vAiCr 

computer systems that were looking at arrhythmias respec- 
tively report finding 71?S or PVG*s,^ 9!^^ of PVC*s,^^ and 
60^ of PVC*s and AFC'b*^^ It is also important to minimise 
the number of false positives In order to preserve the use- 
fulness of the analysis, but this error Is not as critical. 

The second requirement is that the program be fast. 

Sixty times real-time is a useful speed, and a goal that has been 
talked about. Because of the peculiarities of PVC*s and the 
nature of the studies that involve PVC detection, huge 
amounts of data are generated. As mentioned earlier, there 
are likely to be very large numbers of 2k hour tapes that 
will need to be analysed. Real-time analysis would be 
uneconomical in equipment tied up, personnel required, and 
time needed before output data was available. Because a 
complete analysis is not needed, only a PVC detection, there 
should be considerable room for speeding up the process. 

One computer arrhythmia detector has mn as fast as 4 times 

12 

real-time on a PDP-4 and 15 times on a FDF-7, 

Another requirement for a successful Implementation Is 
that It be dependable and easy to use. The personnel using 
the system are likely to be non-technical and/or non- 
medloal. Therefore the set-up, the switching from one tape 
to another, the output, and the record-keeping must all be 
designed to be handled in a simple routine. Any parameters 
that need adjusting. Including speed changes, must be easily 
accessible and held to minimum. The program itself should 
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handle any abnormalities that are at all likely to occur, and 
when it Is unable to handle some problem It should provide 
an indication to the operator of the difficulty. 

The final general requirement for a good system is 
that it be flexible. Not only must it be able to run at 
different speeds, but it must be adjustable to handle widely 
varying, including noisy, data. As more data is gathered 
and analyzed, the parameters should be adjustable for opti- 
mum discrimination, A new algorithm for a module such as the 
detector block should be easily switched in, and the overall 
program should be expandable to include features such as 
statistical anlysis. So, a successful system must be a fast, 
accurate, flexible, easy-to-use method of detecting premature 
ventricular contractions. 

Processing Required 

There are six main blocks that nearly any solution 
would have to encompass. In order, they are input, data 
reduction, detection, measurement, decision, and output. 

Input 

In order to store the BCG in the computer, the analog 
signal must be sampled and converted to digital form. A 
variety of methods and sampling rates, and a variety of 
reasons to support them, are available in the literature. 

The American Heart Association says that 500 samples per 

8econd(sps) and a 9-bit digital output are needed to fully 

6 7 

analyze an ECG, Some groups agree with this rate;*' while 

others feel 200 sps is sufficient,^* and others feel that 

even 120 ' or 125 sps is all that is needed. Not as much 
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detail is required ror a rvG ae-cecxor as xor a xuii exami-» 
nation of the complete ECG, so a lower sampling is probably 
pennissable. 

Data Reduction 

Any sampling rate comppxable to Ihose indicated above 
provides too many points to be handled at very high speeds. 

For example, 200 sps real-time extrapolated to 60 times 
real-time would req\tlre 12,000 individual points per second. 
Eighty usee, is far from enough time to average on each 
piece of information. Therefore, theamount of data must be 
reduced, but without losing important information in the 
waveform. Between QRS's, and that is the vast majority of 
the time, the sign.al is quite flat and with little informa- 
tion value. A good reduction would preserve a QRS in its 
entirety and describe the baseline with only a few nirntbers. 

One algorithm that does this well is the AZTEC preprocessor 
developed by a group at Washington University. That pro- 
gram approximates an analog signal with a series of hori^i^ontal 
line segments and slopes. As long as the input signal stays 
within a preset aperture of the last output value, the output 
stays constant (horizontal). When the threshold is exceeded, 
a new line segment is started and the process repeated. 
Consecutive short segments, i.e. steep paxrts of the wave, 
ere further approximated by a single sloping line. The size 
of the aperture regulates the amount of data reduction and 
the amount of detail preserved. Low amplitude variations 
are lost, while steep slopes are preserved almost completely. 
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They estimate an average data reduction on £CG*s of lOil 
with little loss of information. A similar prcprooeosor* 
which involves an approximation using only horisontal line 
segments, is called a Zero Order Interpolator and detailed 
in another source. This is faster, preserves steep 
segments more faithfully, and has a larger number of output 
segments than AZTEC. The data storage is also simpler. 

Using this interpolated data, an analysis Should be much 
swifter. 

Detection 

First, it must be decided what one is looking for and 
then it must be decided how best to look for it. For a FVC 
analysis* the QBS complexes provide sufficient information. 

The problem at this stage is then to find all QBS's, normal 
and abnormal* and reject everything else. A variety of 
technlqtues have been employed, but they all center on the 
two main distinguishing features of the QRS, its large 
anplltude and large slopes. Some find the beginning of the 
QRS simply by finding a slope of sufficient steepness, ’ ** ** 
another requires a minimum slope to exist for a minimum 
length of time, '' still another requires two slopes that 
eolXootlvely exceed both critical positive and negative slopes 
within 150 msec, to define a QRS,^^ while others make their 
dateotion a combination of sufficient slope and sufficient 
magnitude. In addition to finding all true events, such 
noise problems as steep spikes and baseline shifts should ' 
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should lead to a minimum of false detections. This stage 
Is intimately related to tne next, that of measurement, 
tdeaeurement 

Once the QRS complexes have been found, their critical 

parameters must be calculated for use in deciding if it is 

a PVC. P-waves are not used since they are difficult to 

measure (amplitude typically only twice that of noise ) 

and they are of secondary importance in finding PVC*s. 

Measurements are usually made of the magnitude and width of 

the QHS, the baseline, and the rate or interval between 

complexes. Within this agreement, there is once more a 

wide choice of methods. Some investigators determine onset 

and end of the complex by looking for a flat stretch of a 

certain length-^* or a stretch which has a slope less than 

some cutoff^^ while others invoke procedures utilizing 

12 20 

first and second derivatives. * The magnitude measurement 

may be either a signed number in relation to the baseline 

or a peak-to-peak amplitude. Interval is self<»explanatory, 

but can be chosen to be QQ, RR, or some other regular 

interval. Baseline ^ rot used in deciding on abnormality, 

is used as an Intermediate measurement to provide a local 

reference zero Instead of sn arbitrary constants It can 

be measured as the longest flat stretch, as the magnitude 

at some particular temporal location in reference to the 

QRS, or as the most common value in a histogram plot of 

20 

the total waveform. This is a very critical part of the 
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program* An error in measurement can easily yield values 
that would classify a normal as abnormal , or vice versa* 

Low caid high amplitude and low( baseline shift) and high 
frequency (spikes) noises all contribute their vagaries 
to the problem* peak-to-peak magnitude Is an easy measure- 
ment, but onset and end of the OKS and a good baseline are 
* often obscured because there may be no significant stretches 
of flat wave at the right places* 

Decision 

The decision module is in a sense the heart of the 
analysis, where all the other data is gathered and analysed* 
This project was defined to search only for PVC*s, so this 
algorithm may be comparatively simple. However, in addition 
to deciding wiiat is a PVC, families of r/C*s having the 
same focus and, therefore, siHdlar shape can be found and 
runs or other patterns of FVC*s can be discovered. Features 
of FVC.*s which could be used to distinguish them include 
prematurity, compensatory pause, width, height, bizarre 
shape of QRS, and polarity. The two most widely used are 
width, which one source^ says is the most important number, 
and prematurity. Definitions of prematurity range through 
60j5, 70 ^, 80^, 85?S, and 90^ of normal interval. The 

normal Is generally a running average of the normal intervals. 
This is essential since not only does normal vary greatly . 
from patient to patient, but also greatly for one patient 
over a period of time as long as 24 hours. Some of the 
other features mentioned are also used by some researchers. 
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In general » the parameters have been considered Independently, 
but it is also feasible to accept a certain degree of slrnul- 
taneous variation In each of a number of parameters to 
constitute a basis for abnormality, normally, once PVC*s 
have been found, they would be flagged in some way and per- 
haps processed further in some statistical survey. 

Output 

There Is no particular output which la best in general, 
instead it depends on the specific application. The only 
necessity Is that it provide a means of accessing the 
- abnormals and a coding of why they were considered abnormal. 
The easiest output is to simply supply a running string of 
the parameter measurements and the flags of PVC's. Unless 
delay of some sort is utilized the parameters and flags will 
necessarily appear at the output later than the raw signal. 
This presents no difficulty if It is feasible, in the follow- 
up stage, to play the output tape baclcwards and hence see 
the flag before the FVC. At high speeds it would be difficult 
to provide hard copy output of abnormals in a single pass, 
simply because a strip-chart recorder cannot move that fast. 

It would presumably be not too difficult to build a piece 
of hardware for the second stage, which accepts a magnetic 
tape recording of the £CG, coded parameters, and flags. A 
typical option would be to play this tape backwards, turn 
on a strip chart recorder when a flag is detected, output 
that FVC and as much surrounding context as desired, turn 
the strip chart off, and continue until another PVC is 
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fotmd. If all that is required is a statistical sununary 
or histogram of abnormal events, that can of course be 
provided with only one pass in whatever format the user de- 
sires, Another idea which has been suggested is to make 
use of the digital formatting, and store PVC*s and data on 
digital storage of some kind. A playback of this could be 
done much more flexibly and without a tape recorder. 

The Implementation 

A system which performs the tasks described earlier 
has been established and is working. The actual development 
has been done on a Data General Nova computer. It is 
intended that the present program be transferred without 
change to a Nova 800 in order to utilize better the speed 
potiential of the system. Most instructions on the Nova 
require 5 'to 6 usee, for execution, so that an immediate gain 
in speed of a factor of 4 or 5 could be expected by using 
the faster Nova 800, The Nova is a four-accumulator, 16- 
bit machine. The entire program has been written in assembly 
language, and it currently contains about 500 executable 
Instructions. Ai^sembly language was used because of the 
critical speed rviquirements of the program. A flow chart 
showing the entir^i program in rough blocks is in figure 4. 

The program was developed under a disk operating system, 
permanent storage is on Lino tapes. Figure 5 illustrates the 
toll implementation. Communication between the computer 
and the outside world is done through an Analog Input and 
an Analog Output, which are parts of an Analog Chassis that 
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Figure 4 

Flowchart of master algoritiun 
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was designed to be used In a clinical compuxaxion system. 

They provide j respectively, the buffered analog/digit 2 il 
and digital/su^oE converters. In the same Chassis is a 
High-Resolution Cloch which is program adjustable to pro- 
vide accurate sampling of the signal for any selection of 
processing rate. Communication between user and program is 
via teletype under DOS control. A tape recorder is needed 
for the input signal and the same or different tape recorder 
and/or strip chart recorder are needed to receive the output. 
It is highly recommended that a display scope be provided 
to ascertain that gain and offset on the input are appro- 
priate, and special signals are provided for an X-Y plotter 
if that is desired. 

Input 

The £CG signal from a tape recorder Is delivered directly 
to the Analog Input(AXN) mentioned above. The full input 
range of the AIR is adjustable from ^0.2 volts to 
+<1- volts, and a DC offset of volts provided at the front 
panel. When triggered by the High-Resolution Clock(HRC) 
described Immediately below, a 10-blt analo&^digital 
eonver£bn is made. The period of the HRC can be adjusted 
by the computer in steps of 1 nsec. A 200 sp:? rate has 
been used for real-time analyslfi. so that the HRC provides 
a conversion command to the AIR every 5000 Aa the 

processing rate is increased the HRC period Is decreased 
accordingly, by changing the value of a single variable in 
the program. Every time the AIN completes an Input, an 
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data is stored in a circular buffer that Is long enough 
to store a si^^i4A njbl4> fill the display^ screen* At 

this tine, the data reduction algorithm also considers these 
points t as described following. 

Data Reduction 

The method of data reduction is the sero order interpo- 
lator (ZOX) described earlier* Basically, a horisontal 
line segment approximation is made to the input signal. 
Whenever the signal differs from the most recent output 
by more than a prescribed aperture, a new output segment 
is begun* Alternately, if the output remains constant for 
greater than a prescribed lexigth of time, a new line is 
started in order to provide a reasonably regular data 
stream. A flowchart of thi^ structure is in figure 6* The 
line segments are stored as word pairs, representing the 
amplitude and length of the ;’,egments, In a circular buffer 
which is again sufficiently long to provide a full-screen 
display and to temporarily buffer heart-rate changes and 
lags in the analysis* The sise of the aperture may greatly 
affect the shape of an interpolated signal. A larger 
window will ignore low amplitude deviations, but preserve 
large ones, in addition to providing a larger data reduction* 
Figures 7 through 16 illustrate the properties of the ZOl 
and the effect of changing the aperture* Also illustrated 
is the effect of changing the sampling rate from 200 sps 
to 120 sps* Three PVC*s are included for reference. 
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Flowchart of Zero Order Interpolator 
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Figures 7 through 16 were all produced on an X-Y 
plotter, using a feature of this program. The top 
line is the sampled and digitised signal and the 
bottom line is the output of the zero order Inter-^ 
polator. There are 512 points per line, Bx? that, 
at 200 samples per second, this represents 2.56 
seconds. This is exactly the same picture as seen 
by the operator on the display scope. Unless spec- 
ified otherwise, sampling rate was 200 times per 
second and the aperture was 32 units out of a possible 
1024. The normal QRS was adjusted to occupy approx- 
imately 512 units, so that the aperture represents 
about 6^ of a normal QRS height. 


Figure 7i Interpolation with default parameters. 

Note that the QRS is virtually unchanged by the 
approximation, while the baseline and smaller waves 
have the number of segments representing them greatly 
reduced. It can be easily verified that the number 
of line segments Is about one tenth the number of 
original bits. 
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Figure 8t Interpolation with aperture twice normal. 
This is the same waveform as figure 7, with the 
apertiire now set at about 12^ of the QRS height. 
The P-wave has been lost, the T-wave is still 
easily recognizable, and, except for the very 
onset, the QRS complex is still nearly the same 
as the sampled complex. 







Figure 9i Interpolation with aperture 4 times normal. 
This is again the same waveform as in figiurea 7 
and 8, with the aperture now set to 128 uni'^s. 

The T-wave has lost its morphology, although 
its looation is well marked. The QRS is still 
very well represented. The data reduction is now 
about 15»1. 
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Plgiire 10 1 Interpolation with normal aperture , at 
120 samples per second. This is the same waveform 
and aperture as figure 7, with 3/5 the number of 
sampling points per second. The width of the 
picture is now 4,27 seconds. In the upper trace, 
the amount of detail is noticeably less than at 
200 samples per second, but probably quite sufficient 
for analysis. Again, the interpolation at this 
aperture provides an excellent reproduction of the 
QRS. 
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Figure ill Interpolation with aperture twice normal, 
at 120 samples per second. This corresponds with 
figure 8 in its data reduction and feature repre- 
sentation. The P-wave is gone, the T-wave 
detectable, and the QRS totally reproduced with 
the exception of the tip of the Q-wave. 
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Figure IZt Interpolation of a noisier recording, 
default parameters. This tracing was done from 
data originally gathered by American Optical. 

The large amount of high frequency noise is 
inadequately represented by the sampled signal. 
The noise is reduced somewhat by interpolation, 
but there remain a large number of short segments 
in the baseline. These present a problem as far 
as degree of data reduction is concerned, and as 
far as recognizing legitimate steep slopes is 
concerned. 



Figure 131 Interpolation of a noisier recording, 
with aperture twice normal. This is the same 
waveform as figure 12, with the aperture at twice 
the coarseness. The baseline Is now very reasonable, 
but note that many of the segments in the QRS are 
now two units long. This is not necessarily a 
problem, but must be taken into consideration by 
the detection algorithm. 



Figiire 14* Interpolation of a PVC, with aperture 
twice normal. This display is from the same 
recording as figures 12 and 13 . Note that both 
the PVC and the raormals stand out well from the 
filtered haseline, but that again care must be 
talcen to ensure that segments within the QR3 of 
lengths greater than 1 unit are properly inter- 
prettcd, 



Figure 15t Interpolation ot a second FVC, default 
parameters. This waveform is from a tape recorded 
at Peter Bent Brigham Hospital, The very clean 
signal depicts well the features of thes algorithm, 
The negative portion of the PVC is flattened due 
to saturation of the analog/dlgltal converter. 



Figure l6t Interpolation of a third PVC, default 
parameters. This is a different FVC from the 
same PBEK tape. Again, the necessary features 
are all well-defined, but the detection 
algorithm must be carefully designed to abstract 
the proper conclclusions about onset and end of 
the waveforms. 


50 




* ■«- ■ • *. 




•'V 




% 




OfoOmiAi; PA6B IB 
MEtX»{inAURi 


51 


Currently, the default value of this aperture is about 6^ 
of a normal QRS peak-to-peak amplitude. As can be seen, 
the QRS is preserved nearly exactly, all significant wave- 
forms are still easily recognizable, and the nximber of line 
segments is approximately a factor of 10 lower than the 
number of points input. The sampled signal and the output 
of the ZOI are displayed on the display scope. The figures 
included are all drawn directly on the X-Y plotter, using 
the option described in the section on how to use the 
routine, and are the same pictures that appear on the display 
scope. 

Detection 

This part of the program attempts to find all normal 
and abnormal QRS complexes, and reject everything else. 

Its basic premise is that QRS*s are comprised primarily of 
large slopes, most large slopes are in the QRS, and the QRS 
is larger than the rest of the EGG. Initially, the program 
waits until a long flat section has passed, to ensure that 
its detection begins between complexes. When a segment of 
oiritlcal slope is found, it is labelled as a possible QRS 
onset. As long as the signal remains steep, it is considered 
part of a possiole QRS, A certain length of time after 
onset, the peak-to-peak amplitude so far observed is compared 
with a fraction of the running average of the total peak-to- 
peak amplitudes. If it is smaller, the event is rejected 
as being too small for a QRS, and a new search is begun. 

When a flat stretch of specified length is detected, it is 
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considered a possible QRS end. If it is further from the 
average baseline than a certain fraction of the average 
peak-to»peak amplitude, that flat stretch is considered to 
be an elevated part of a PVC, and the search for a QRS end 
continues. Otherwise, it is decided that the full event has 
been detected. If the event is narrower than a specified 
cutoff, it is termed a noise spike, otherwise it is a valid 
QRS complex. Also, all QHS's are automatically terminated 
after a certain length of time (such as 200 msec.) if they 
have not already been by a flat segment. The flow chart 
in figure 17 illustrates this algorithm. After the end 
of the QRS, a certaiin length of time(150 ms. default) is ignored 
so that no large, steep T-waves will be mistaken for QRS 
complexes. 

Reasonable baseline shifts, small amplitude waves and 
noise, and narrow noise are all successfully rejected by this 
algorithm. One critical parameter must be adjusted to dif- 
ferentiate between thin QHS's and noise spikes. On record- 
ings with no more noise than can noinaally be expected in a 
ellnloal environment, this module does an excellent job of 
finding all QRS's and ignoring the other events. There is 
a large room for modification and improvement to upgrade 
the detection under more severe environments. The onset and 
end are not always found precisely on difficult samples, 
but this is more of a measuremnent problem, as explained in 
the next section. The critical slopes, the critical flat 
sections, the times where different parameters are observed, 
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and the amplitude cutoffs for event decision are all single 
parameters that can be adjusted with experience and for 
data that is unusual. The enclosed sample strip charts 
Illustrate some of the successes and failures of this routine. 

Measv*remant 

For each QRS that is detected, certain parameters have 
to be calc\ilated. Their accuracy is heavily dependent on 
the success of the detection module. Width and magnitude 
of the QRS, Q4 Interval between QRS*s, and baseline amplitude 
are measured. In order to make time calculations, an aux- 
iliary variable called SegmentClock has been created. As 
each line segment is read form the circular buffer its length 
is added to SegQlk, to unlciuely (modulo 2 *5 msec.«B5 min. 
at real-time) specify its temporal location. The width is 
simply elapsed time between what was detected as onset and 
end of the QRS, QQ interval is between consecutive onsets, 
and magnitude is determined by a running window which expands 
as each new maximum or minimum within the QRS is found. 
Baseline for a particular complex Is defined as the amplitude 
of the last flat segment prior to the complex. If, as in a 
very large PVG,the input signal saturates the analog/dlgltal 
converter, the magnitude is calculated as if the signal had 
barely touched the extremum. The measurement routine is 
consistent, as observed on several identical analyses of the 
same data. The maximum width measured is determined by the 
detection routine, now set at 200 msec. The program halts 
if it finds a QQ interval less than 240 msee. or greater 
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than 10 seconds. The variables that affect measurement are 
those mentioned above in detection, primarily for making 
accurate decisions on where QRS onset and end are. The main 
cases of Inaccuracy are steep baseline shifts and noise 
which camoflages onset and end. Coarsening the interpolation 
process does help with low amplitude, high frequency noise. 

Shortening the flat stretch indicating an ST segment allows 
accurate detection of end where it would otherwise be 
missed in noise, but it reciprocally introduces a hazard 
of terminating too early. Permitting a flat segment far 
from the baseline to be called a OHS end instead of a flat 
segment of a PVC filters out baseline shift better but 
increases the probability of mls-measuring PVC parameters. 

Decision 

Presently, only binary decisions are made on PVC versus 
normal QHS. All data is available to classify PVC*s into 
families having similar morphology if that is so desired. 

To find abnormal beats, the measured parametas are compared 
to running averages deteitnined from the normal beats of 
that patient. The relative v/eighting of the averages is 
adjustable, presently the most recent parameter is given 
1/3 or 1/7 the influence of the previous average. A PVC 
is found when the current width is significantly larger 
than the running average, or the QQ interval is significantly 
shorter. The magnitude is not considered In this algorithm. 

The actual cutoffs are easily adjustable, with defaults at 
175^ of average width and 15/15 of average QQ. These are 
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yields a better foundation and as a function of the indi- 
vidual patient *s i^.uu. For example, a person whose normal 
heart rata varies considerably with breathing v;ould want a 
less sensitive QQ cutoff. It is esential that running 
averages be used as the standard, since normals vary widely 
between patients and over time in the same patient. The 
.first several complexes in each new tape are automatically 
considered normal, in order that their parameters may be 
incorporated into an initial running average for that tape. 

A more complex decision tree could be easily implemented. 
Maginltude of the QRS and information on compensatory pause 
are accessible, and a decision that considered the depen- 
dency of variables would also be feasible. This module 
is the main area where tradeoff between false negatives and 
false positives is determined. Tighter bounds will let fewer 
abnormals by, and also classify more normals as abnormals. 

The current fairly loose bound on the width is due to the 
fact that measurements of width are not highly precise. So 
far, the decision algorithm has proved accurate, with blame 
for error resting primarily on detection and measurement. 

Output 

The output consists of a single channel where the 
amplitude, width, and QQ interval appear coded as three 
pnlsea whase height is proportional to the value. Another 
pulse, opposite in polarity to the first three, appears when 
there is a FVC. The distance between sets of pulses Is 
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equal to tne uvt inxerval minus the total wicUh ox the pulses 
so that the output remains a constant time delay behind the 
input ECG. This delay is adjustable. The output parameters 
are stored in a circulsir buffer to allow for variable delay 
and a certain amount of variation in speed of analysis* It 
is expected that this data channel would either be recorded 
on a separate track of the ECG tape, or the ECG and the data 
would be recorded on 2 channels of an output magnetic tape. 

If all that is desired as an output is a set of internally 
generated statistics on PVC’s, there is no need to record 
the data channel. Otherwise, It Is intended that the coded 
channel would be usable for further analysis and to enable 
the PVC*s to be accessed very easily, such as for ebtaining 
strip chart recordings of only the abnormal beats. More 
Ideas on this are presented in the paragraph below on future 
options* Figures 18 through 27 illustrate many features and 
problems of the system as recorded from clinical data* 

Further explanations are included with the recordings. 

How to Use 

The use of the total program is very simple* First, 
the sampling rate must be changed If is to be other than 
real-time (or whatever other permanent rate has been installed). 
If any other parameters are to be changed from default values, 
that should also be done at that time. The four sense 
switches should be off, for the display mode. The tape 
recorder is attached to the Analog Input, and run at the 
appropriate speed. A dual trace now appears on the display 
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Figures 18 through 2? were produced on a strip chart 
recorder with tape speed 25 mm/second unless noted 
otherwise, and sensitivity 100 mv/di vision on both 
channels in all cases. The EGG waveform is the raw 
signal from the tape recorder, while the lower 
channel is the output provided by the computer via 
the Analog Output, The height of the three pulses 
are proportional, in order;*, to the peah-to-peak 
amplitude of the QRS complex, the width of the QRS, 
and the QQ interval. A fourth pulse, of a constant 
negative value, is present in the case that a QHS 
is declared to be a PVC, The data channel is delayed 
behind the QRS by an adjustable amount, in these 
cases about 4 seconds (the full-width vertical red 
lines are 4 seconds apart) , 

Figure I8i A detected PVC. The large PVC at the left 
has been detected, its flag is to right of center. 
Note the similarity of the coded parameters for the 
other CiRS*s, as contrasted to the large amplitude 
and width and short QQ interval indicated for the 
PVC, The magnitude is not considered in the current 
decision algorithm, but either the width or the 
Interval would have classified this as a PVC. The 
following set of parameters indicates a large QQ 
interval, the compensatory pause. The delay between 
PVC and flag is the appropriate 4 seconds. This 
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was recorded from a PBBH tape. 
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Figure 19 x Another detected PYC. As in figure 
18, the PVC was recorded as being large in 
magnitude and width and short in QQ interval. 
The steep P-waves posed no problems. This is 
from the same PBBH tape, and the QQ cutoff was 
set at 7/8 of average. 
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Figure 20i Two more detected PVC*s. These two strips 
form a continuous segment of tape. At the beginning 
of each strip is a PVC, plus an abnormal non-PVC 
near the end of the first strip. The two PVC's 
were flagged as such, theother was not. Note the 
unsteady baseline. The abnormal non-PVC passed 
because it was not sufficiently premature, and its 
measured width was about normal. Again, the PBBH 
tape and QQ cutoff of 7/8 average were used. 
















Figure 21 I A longer recording. This is a contin- 
uous recording made at 5 mm/second to allow a 
visualization of a longer pattern of detection. 
The only 2 PVC's were detected easily, and there 
were no false positives. It can be seen that 
the measurements are reasonable, and do follow 
slight variations such as magnitude from beat to 
beat. This is from the PBBH tape. 
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Figure 22i PVC*s from another patient. This 

recording is from the American Optical tape. The 
QQ cutoff is 7/8 of average interval, the inter- 
polator aperture is twice normal, and the critical 
length to determine a QRS slope has been raised 
from 1 to 2 (in conjunction with the coarser Inter- 
polator) . This continuous recording contains 2 
PYC*s plus fair amounts of high frequency noise. 

In this case the peak-to-peak amplitude of the PVC*s 
are less than for the normals, the widths are of 
the same order, but the QQ intervals are well within 
the cutoff of 7/8 of the running average. The noise 
has not significantly affected the measuremnets. 
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Figure 23i A direct recording at 8 times real-time. 
For this sample, the PBBH tape was run at 8 times 
normal speed (and the sampling rate increased 
8 times) , and the strip chart was speeded up 5 
times to 125 mm/second. The frequency response 
of the strip chart recorder is evident in the 
deterioration of the pulse shapes, hut the algo- 
rithm and output work fine. 
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Figure 24i An indirect recording at 8 times real- 
time. In this example, the same tape was again 
played at 8 times the normal rate. But this time 
the output was recorded on a second magnetic tape 
(as would he expected in normal performances). 
This magnetic tape was played back at \ ox its 
recording speed and the signal fed to the strip 
chart recorder, which was played at 25 mm/seconc\. 
The analysis is the same as in figure 23, but the 
paper recording is enhanced. 
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Figure 2$t A false positive identification. This 
recording is from the American Optical tape. On 
the third QRS, there is a significant- amount of 
steep noise at the S-wave. Due to this the end 
of the QRS was measured later, the complex was 
therefore incorrectly called vdde, and the beat 
flagged as a PVC, Also, immediately prior to the 
fifth QRS is a negative, narrow noise spike. It 
can be seen that the program realized that it was 
a spike and reset in time to identify the Following 
QRS correctly. The aperture is twice normal, the 
length to determine QRS onset is 2 instead of 1, 
and the QQ cutoff is 7/8 of an average interval. 
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Figure 261 A second false positive identification. 

In this case the large amplitude, steep noise was 
assumed to be a QRS complex. The amplitude was 
large enough, and the 2 largest spikes may have 
been observed as a single complex to yield an 
appropriate width. This p«f»ndo—cntnr>i wap 
considered premature and therefore the following 
resLl complex also had. to be assumed premature. 

Both were flagged. The aperture was 1.5 times 
normal, the QQ cutoff was 7/8 of an average interval, 
and the amplitude needed to declare an event a QRS 
complex was half of the average amplitude. This 
is from the American Optical tape. 
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Figure 2?t A third false positive identification, 
and its solution. The EGG sample for both strips 
is exactly the same segment from the PBBH tape. 

On the top strip the QQ cutoff is 7/8 of an 
e.verage interval. On the lower strip this is 
also true, and in addition the fiat length needed 
to assume an ST segment has occurred is 30 msec, 
instead of the default 50 msec. On the top strip, 
the large swing in the baseline was sufficient to 
confuse the identification of QRS end. Therefore 
the QRS was assumed wide and flagged as a PVC. 
This problem was circumvented by allowing the QRS 
to terminate on a shorter flat segment, as can be 
seen in the lov;er tracing. 
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scope. On top is the sajnplecL £.CG, ana on the ootxoin is xne 
output of the zero oraer interpolator. The screen is con- 
tinuously refreshed so that the traces appear to move slowly 
to the left. Figure 2b contains a flowchart of the display 
process. If sense switch u is raised in either the display 
or evaluation mode, the process returns to DOb control. 

Raising sense switch 2 freezes the display, so that a 
detailed look or pnotographs can be made. Lowering the 
switch returns the moving display*. While sense sw^.tch 
2 is up, raising sense sv/itcn ^ will enable the X-Y plotter 
and then halt so that paper can be loaded, etc. When 
the continue switch on the computer is pressed, a sweep is 
made on the X-Y plotter of whatever had been frozen on 
the display scope. As long as switch 3 is up, the same plot 
will be repeated. Lovrering switch 3 returns the program 
to frozen display mode and disables the X-Y plotter. The 
figures included earlier were done this way. 

The primary purpose of displaying the signal as the 
computer sees It is to allow the gain and offset to be 
adjusted to values which maximize information. It is recom- 
mended that normal peak-to-peak magnitude be adjusted to 
occupy about half the screen. When this has been done, sense 
switch 1 is raise and the program begins evaluation and 
output. The display stops, as this is very time consuming. 
Once the program is running, it will do so until a QQ interval 
is less than 240 msec, or more than 10 seconds, until an 
I/O device does not respond properly, until the program is 
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Uiiaule to lvaa*dlG tt^a data rate as indicated "by the data buffer 
being emptied > or imtll it is inanualXy stopped, Restarti3ig is 
easy, using the debugger, and is essential for a new tape because 
new ruiining averages have to be establiched. If the program is run 
under debugger control, all features of that mode are available. 

Theoretical Speeds, Rtuming Times 

At real-time, on the Nova computer, the EGG is sampled 
2Q0 sps, or eve2?y 5 msec. To handle each sample of input, 
about 50 instructions occupying 30O usee, are required, 

Therefore, about 60 msec, out of everj^ second is needed to 
handle the input stream. This would indicate a maximum 
of 16 times real-time, for just inputting the signal, 
executing the ZOI, storing the data, and outputting the 
coded parameters and flags. To estimate the time needed 
for analysis of the wave, a heart rate of 80 beats per minute, 
or oeat per second, is assumed, A generous estimate 
is to allocate l^ segments to the QHS complex and another 12 
segments to represent the interval between complexes. 

Approximately 450 instructions are executed bety/een each 
complex, 720 are used in analysing the QRS, and about 250 
arc needed for handling bookkeeping after each complex. 

Therefore there are about 1420 instructions executed per 
complex and I900 per second, still assuming real-time. This 
represents about 11 msec, out cf each 1 second of data. 

The total time required by the program is hence about 
71 msec, out of each second at real-time. As can be readily 
seen, the time spent on analysis is a small fraction of the 



time spent on 1/0 handling. Therefore the estimate of 
71 msec, is relatively insensitive to heart rate or number 
of PVO*s, hut rather primarily a function of the sampling 
rate. Using these numbers, a program should be able to 
run at times real-time. If this is true, minor changes 
could be made in the instruction coding to increase the 
rate to 16 times real-time without changing the algorithm 
at all, 

A large increase in speed would necessitate a decrease 
in sampling r? ' For 120 sps, only 3/5 ot the amount of 
data is hand! The 14 times real-time would then go to 
about 2i tim< Implementing any verson on the Kova 800 
is expected to yield an automatic increase in speed by a 
factor of 4 or 5 * With this in mind, rates of 60 times 
real-time for the 200 sps model and 100 times for the 120 
sps model are not unreasonable expectations. The 120 sps 
would of course lose some detail, but that is a rate which 
has successfully been used previously, as was referenced 
in an earlier section on background. Figures 10 and 11 
were made with a sampling rate of 120 sps. 

Analysis of Results 

The four criteria of accuracy, speed, ease of use, and 
flexibility were presented as the requirements for a 
successful high-speed evaluation, of ECG's, Accuracy has been 
established to the extent that it works on several tape- 
recordings which include PVG's. Several dozen sections of 
tape which were known to contain PVG*s have been analyzed, 


and all PVC's found, Very few non-PVC*s have "been labelled 
as such, and most of these false identifications can he 
eliminated with a more knowledgeable adjustment of para- 
meters and criteria. I’he rate of false positives is generally 
below X^t sometimes well below. This number is still higher 
than desired, but should be greatly reduced by using the 
full power of the program. A thorough statistical analysis 
of the accuracy under a variety of conditions is critical, 
and should be the next step in furthering this project. 

Based on experience to date, the accuracy of the program 
is good, particularly as regards having caught every PVG. 

The strip chart recordings included in this paper shov/ 
a variety of PVC*s that have been detected, illustrate 
a couple of situations that lead to false positive iden- 
tification and point out hov/ to eliminate some of these 
errors, and attempt to depict the overall accuracy of the 
algorithm. 

The program has been run without any difficulty at 8 
times real-time. Figures 23 and 2lp illustrate this. 

Attempts to run at 16 times r-sal-time start successfully 
but (luickly halt when the demand for output exceeds the 
speed of analysis. This would seem to support the theoret- 
ical estimate of about 14 times real-time given above. 

This speed can be Increased significantly as was described 
in the section on theoretical speeds. Although not 
ouXMtanding, it seems reasonable that this operation can 
be run at 60 times real-time or more on ths Nova 800 , as 
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originally intended. 

Ease of use is certainly inherent in thits system. If 
no parameters are to he adjusted, the total operation consists 
of starting the tape recorder, adjusting the gain and offset 
as monitored on the display scox©, and raising one switch 
to enter evaluation mode. Changing parameters is also very 
simple under debugger control, 

^Dhe entire program is q[uite flexible in regard to 
modularity, changeability of pairameters, and ease of adding 
new options. Entire blocks such as the detection and 
decision schemes may be replaced by others if so desired. 

As explained many times, individual parameters are easily 
changed, and should be to obtain better accuracy as better 
infoimiation is received and on patients with v;idely 
dissimilar EGG'S, She present output provides a basis that 
can be used for a variety of output options, some of which 
are explained belov/ under suggestions for the future. 
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SUGGESTIONS FOR THE FUTURE 


This basic program can be improved in many ways, can 
be modified for similar piirposes, and can have feattiree 
added to increase its power. The first task should be to 
make a clinical, statistical evaluation on a large bulk of 
data. Presumably, the detection and decision blocks in par- 
ticular have ample room for improvement. If that amount of 
increase of speed is desired, the idea of changing the 
sampling rate from 200 to 120 sps should also be tested. 

This algorithm has been specifically programmed for PVC 
detection, but there is no reason that certain other 
arrhythmias could not be studied using a different set of 
parameters in the same basic framev^ork. Even for PVC*s, 
a completely different detection algorithm could be created. 
For instance, the initial detection could be based on 
amplitude instead of slope* 

The output provides for the largest variety of options. 
This has been delibera±,ely left open-ended, since each user 
will probably have different requirements. The program 
could be amended to do a statistical analysis of the input 
signal and provide a teletype output at the end of each 
tape. Or, the normal data channel output could be created 
and t^sed as an input to a second stage hardware system 
that displayed all the PVC’s on strip chart. Or, the program 
could store its output in digital memory v;ithout ever pro- 
viding an analog signal. This idea could be very flexible 
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in iiis implementation, and also very power rui, TJtie system 
could also Tae structured so that the operator could select 
a particular output option as a function of the result of 
the evaluation. 
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GOWCLUSION 


There are currently several applications where it is 
desirable or even necessary to evaluate magnetic tapes of 
ECG*s at many times real-time. In particular, PVC»s are 
of great clinical and research interes*^. Much work has 
been done on computer analysis of ECG’s in general, but 
this application has not been well-investigated. 

An effective solution must have high accuracy, be 
rapid, be flexible,, and be easy to use. The general blocks 
that are required in such a program are input, data reduc- 
tion, detection, measurement, decision, and output. 

Such a system has been implemented and works, on a Nova 
computer. On a clinical data base including over 75 PVC*s, 
every PVC was flagged. The number of normal beats and 
artifacts identified as PVO’s was under Vfo, The tested 
speed is over 8 times real-time and has been calculated to 
be about times on the Nova. Implementation of the 
same program on the Mova 800 will increase the speed by a 
factor of or 5* Flexibility and expandability are 
inherent in the design. 

At this point, a thorough analysis of the effectiveness 
of the system must be done, on a larger data base. In 
addition to determining the actual power of the system, 
this will provide the basis on which to establish better 
values for the program parameters. It Is es^sected that 
this step will greatly reduce the rate of false positives. 


One option that could increase the speed of analysis 
by 50 ^ is to reduce the sampling rate from 200 sps to 120 
sps. The output hloch may he changed or added to in 
order to mold the system for a particular application. 

It is hoped that other changes will also he made without 
affecting the overall flow of the program, to elevate 
this from a v/orhing system to a useful tool. 


88 


REFERENCES 


1, Black, Herbert, "Medical PulseA'^achines Monitor Heart 
Constantly", Boston Globe . c.April 1, 1973* 

2, Bonner, R.E, , "A Computer System for EGG Monitoring", 

IBM Technical Report no, IBM Corp., Advanced 

Systems Development Divis ion, Yorktovm Heights, K.Y, 

• 3* Caceres, Cesar A. and Dreifus, Leonard S,, Clinical 

Electrocardioaya'Dhv and Qonmuters , Academic Press, Inc,, 
1970, rTew York, K,Y, 4? 1 pages. 

A, Cox, J.R.Jr. , et, al, , "Some Data Transformations Useful 
in Electrocardiograph^' , Gomnuters in Biomedioal Research , 
ed, Stacy and Vfaxman, vol.3. Academic Press, Nev/ York, 
pp. 181-206, Quote was from p.l9^* 

5. Cox, J.R.Jr, and Nolle, P.M., "Arrhythmia Monitoring 
Algorithm for Real Time Application" , Proceedings of th e 
Hawaii International Conference on System Sciences^* 
Computers in Biomedicine . ■dp, 120«122, Honolulu, Hav/aii. 
1972. 

6. Cox, J.R.Jr, , et.al,, "Digital Analysis of the EEG, the 
Blood Pressure Wave, and the EGG", Proceedings of the 
IEEE . v.60 no. 10, October, 1972, pp. 11 37 ”11 64, 

7. Cox, J.R.Jr,, et.al,, "A2TEC, A Preprocessing Program for 
Real-Time ECG Rhythm Analysis" , IEEE Transactions on 
Bio-^edical Engineering . v,15 mo, 2, April 1968, pp, 128-9, 

8. Feldman, C.L., et.il., , "Computer Detection of Ventricular 
Ectopic Beats", Computers in Biomedical Research . v.3» 
p,666, 1970 , 

9. Kohn, M.S., "Characterization of Ectopic Heartbeats", 

SM Thesis, Electrical Engineering, MIT, 1968 , 

10, Kumik, P.B., "Computer Aided Analysis of Electrocardio- 
grams", Tech Engineering News . v«54 no.7» December 1972, 
pp. 7 - 16 , 

11, Lown, Bernard, "Intensive Heart Care", Scientific American . 
V , 219 no.l, July 1968 , pp, 19 - 27 , 

12, Lynch, J,T. "Arrhythmia Detection and Classification 
Using a Digital Computer on Electrocardiographic Data" 

SM Thesis, Electrical Engineering, MIT, June 1965 * 


89 


13. Naylor, W.S., et.al., "A Reliable ECG R-Wave Detector 
Applicable to_ Normal and Abnormal Rhythms" . Pi goat of 
the Seventh International Confartirise on y;elicL>l .^'id 
Biological Engineering . p,306, Stockholm, Sweden, 196?, 

Nolle, P.M., Cox,J.R.Jr. , "How to Keep Up With the ECG". 
Digest of the 1969 IEEE Comiouter Grouo Conference , 
pp, 182-135* 

15* Nolle, F.M., "Arrhythmia Monitor for Cardiac Patients", 
24th Arurmal Conference on Engineering in Medicine ^d 
Biology . Las Vegas, Nevada, p.l50, 1971* 

16, Pryor, T.A., et.al., "Electrocardiographic Interpretation 
by Computer", Comnuters and Biomedical Research . v,2, 
pp. 537-5^8, 1969* 

17* Schwann,H,P. .editor, Biological Engineering . Gesselov/itz , 
D.B* and Schmitt, O.H., "Electrocardiography", pp*333- 
391, McGraw-Hill, 1969. 

18, Walters, J.B.Jr., "Portable Clinical or Laboratoi-y 
Computing System", Quarterly Progress Rsuort . No.106, 

July 15, 1972, MIT Research Laboratory of Electronics, 
pp, 172-184. 

19, Walters, J.B.Jr*, "Signal Preprocessing as an Aid to 
On-Line EKG Analysis", SB Thesis, Electrical Engineering, 
MIT, 1971. 

20, Vtotak, Josef, Comnuters in Electrocardiogranhy . Charles 
G, Thomas, Springfield, Illinois, 1970, Quote was from 
p.vii. 


